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ABSTRACT 

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg^ 
of imaging with photometry in 21 bands extending from the far-UV to the far-IR. 
These data complement our spectroscopic campaign of over 300k galaxies, and are 
compiled from observations with a variety of facilities including: GALEX, SDSS, 
VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by 
VST and scheduled for observations by ASKAP. These data are processed to a com¬ 
mon astrometric solution, from which photometry is derived for ^ 221,373 galaxies 
with r < 19.8 mag. Online tools are provided to access and download data cutouts, or 
the full mosaics of the GAMA regions in each band. 

We focus, in particular, on the reduction and analysis of the VISTA VIKING 
data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the 
data and examining its integrity. Having derived the 21-band photometric catalogue we 
proceed to fit the data using the energy balance code MAGPHYS. These measurements 
are then used to obtain the first fully empirical measurement of the 0.1-500/im energy 
output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) 
across three time-intervals (0.3-1.1 Gyr, 1.1—1.8 Gyr and 1.8—2.4 Gyr), we find that 
the Universe is currently generating (1.5 ± 0.3) x 10^^ hyo W Mpc“^, down from 
(2.5 ±0.2) X 10^^ hyo W Mpc“^ 2.3 Gyr ago. More importantly, we identify significant 
and smooth evolution in the integrated photon escape fraction at all wavelengths, with 
the UV escape fraction increasing from 27(18)% at ^ = 0.18 in NUV(EUV) to 34(23)% 
at z = 0.06. The GAMA PDR can be found at: http://gama-psi.icrar.org/ 

Key words: galaxies: general — galaxies: photometry — stronomical 

databases: miscellaneous — galaxies: evolution — cosmology: observations — galax¬ 
ies: individual 


1 INTRODUCTION 

Galaxies are complex systems. At the simplest level ionised 
gas cools within a dark matter halo (White V Rees 1978), 
condensing in the densest environments to molecular hydro¬ 
gen (Shu, Adams V Lizano 1987) which may become self- 
gravitating and lead to the formation of a stellar population 
(Bate, Bonnell V Bromm 2003). The stars replenish the 
interstellar medium through supernovae, winds, and other 
mass-loss processes (Tinsley 1980; Schoenberner 1983) lead¬ 
ing to metal enrichment, dust formation, and the heating 
of the interstellar medium through shocks and other turbu¬ 
lent processes (McKee V Ostriker 2007, see also Fontanot et 
al. 2006). 

The dust attenuates (through absorption and scat¬ 
tering) a significant portion of the starlight (Calzetti et 
al. 2000), up to 90% depending on inclination for disc sys¬ 
tems (see Driver et al. 2007) and the internal dust geometry 
and composition. The absorbed fraction of the UV/optical 
light (highly dependent on morphology but typically 30 per¬ 
cent for local Universe disk galaxies) is re-radiated at far-IR 
wavelengths (Popescu V Tuffs 2002; Tuffs et al. 2004; Driver 
et al. 2008). Throughout this process gas is being drawn into 
the galaxy from the intergalactic medium (Keres et al. 2005), 
outflows driven by supernova expel material (Veilleux, Cecil 
V Bland-Hawthorn 2005), and tidal interactions with neigh¬ 
bouring dark matter halos may lead to further mass-loss 
(Toomre V Toomre 1972), or mergers (Lacey V Cole 1993), 
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as well as driving gas to the core leading to re-ignition of 
the central super-massive black hole (Hopkins et al. 2006). 
In short, galaxy evolution is governed by a very wide range 
of complex processes that give rise to multiple energy pro¬ 
duction and recycling pathways traced from X-ray to radio 
wavelengths. 

Traditionally galaxy surveys have been predominantly 
single facility campaigns (e.g., the SuperCOSMOS Sky Sur¬ 
vey and other Digitised Plate Surveys, Hambly et al. 2001; 
SDSS, York et al. 2000; 2MASS, Skrutskie et al. 2006; 
IRAS, Soifer, Neugenbauer k Houck. 1987; FIRST, White 
et al. 1997; HIPASS, Barnes et al. 2001) and as a result 
only capable of exploring a fairly narrow wavelength range. 
Therefore they often only probe one constituent of this 
process, e.g., radio surveys which sample the neutral gas 
content (Barnes et al. 2001), optical campaigns sampling 
the stellar population (York et al. 2000), and far-IR cam¬ 
paigns sampling the dust emission (Soifer et al. 1987). While 
panchromatic datasets of relatively modest size have been 
constructed (e.g., the Spitzer Infrared Nearby Galaxy Sur¬ 
vey, Kennicutt et al. 2003), they are generally too small 
to allow a full exposition of, for example, environment and 
stellar-mass dependencies, or subdividing samples to man¬ 
age co-dependencies. 

Part of the problem in assembling a comprehensive 
panchromatic catalogue is the range of facilities required, 
which in many cases are mismatched in sensitivities and 
resolutions. There are also significant logistical issues: the 
physics underpinning the energy processes at each wave¬ 
length are often very different; the distinct data-streams of¬ 
ten have very different wavelength dependent issues requir¬ 
ing a broad range of specialist skills, and the lack of coop¬ 
erative global structures to coordinate observations across a 
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suite of facilities which cross international borders. Sampling 
the full energy range therefore requires cooperation and col¬ 
laboration across a number of subject areas, the cooperation 
of time-allocation committees, extensive resources to man¬ 
age the many data-flows in an optimal way, new techniques 
to combine the data in a robust manner, and an open skies 
policy towards final data-products by national and interna¬ 
tional observatories. 

Progress in this area has mainly been driven by tech¬ 
nological advancements, coupled with large collaborative ef¬ 
forts, and predominantly in two ways: (1) the construction of 
increasing samples of well-selected nearby galaxies, often on 
an object-by-object basis across the wavelength range (e.g., 
the Atlas of SEDs presented by Brown et al. 2014 and the 
S^G collaboration which now samples over 2000 galaxies, 
see Sheth et al. 2010 and Munoz-Mateos et al. 2015); or (2) 
the concerted follow-up of the deep fields observed by the 
Hubble Space Telescope (e.g., the HST GOODs, Giavalisco 
et al. 2004; HST GOSMOS, Scoville et al. 2007; and HST 
GANDLES, Grogin et al. 2011 and Koekemoer et al. 2011, 
in particular). In the former the sample sizes are modest 
(~100s—1000s of objects), in the latter the galaxies sam¬ 
pled are predominantly at very early epochs (i.e., z > 1). 
In short no highly complete panchromatic catalogue of the 
nearby galaxy population exists, suitable for comprehensive 
statistical analysis, while also covering the full energy range. 

The Galaxy And Mass Assembly survey (GAMA; 
Driver et al. 2009, 2011; Baldry et al. 2010) is an attempt to 
provide a comprehensive spectroscopic survey (Robotham 
et al. 2010; Hopkins et al. 2013; Liske et al. 2015) com¬ 
bined with comprehensive panchromatic imaging from the 
far-UV to far-IR and eventually radio. Results to date are 
based mostly on the spectroscopic campaign combined with 
the optical imaging to explore structure on kpc to Mpc 
scales, in particular the GAMA group catalogue (Robotham 
et al. 2011), the filament catalogue (Alpaslan et al. 2014), 
and structural studies of galaxy populations (e.g., Kelvin et 
al. 2014). 

Here we introduce the panchromatic imaging which has 
been acquired, by us or other teams, over the past five years 
from a variety of ground and space-based facilities. These 
surveys collectively provide near-complete sampling of the 
UV to far-IR wavelength range, through 21 broad-band fil¬ 
ters spanning from 0.15—500/im. The filters represented are: 
EUV, NUV, ugriz, ZYJHKs, Wl, W2, W3, W4, lOO/xm, 
160/im, 250/im, 350/im, and bOOgm. The contributing sur¬ 
veys in order of increasing wavelength are: the GALEX 
Medium Imaging Survey (Martin et al. 2005) plus a dedi¬ 
cated campaign (led by RJT), the Sloan Digital Sky Survey 
Data Release 7 (Abazajian et al. 2009), the VST Kilo-degree 
Survey (VST KiDS; de Jong et al. 2013); the Vista Kilo- 
degree INfrared Galaxy survey (VIKING; see description of 
the ESO Public Surveys in Edge et al. 2013), the Wide- 
held Infrared Survey Explorer (WISE; Wright et al. 2010), 
and the Herschel Astrophysical Terahertz Large Area Survey 
(Herschel-ATLAS; Eales et al. 2010). All of these facilities 
have uniformly surveyed the four largesfj^ GAMA regions 


^ GAMA’s hfth region, G02, covers 20 sq deg and overlaps with 
one of the deep XXM XXL helds, see Liske et al. (2015) for further 
details. 


referred to as GOO, G12, G15 and G23 (with only the latter 
held not covered by SDSS). In the future the GAMA re¬ 
gions will be surveyed at radio wavelengths by ASKAP (as 
part of the WALLABY or DINGO surveys) and at X-ray 
wavelengths by eROSITA. 

Gombined, the four prime GAMA regions cover 
230 deg^ and have uniform spectroscopic coverage to 
^Petro < 19.8 mag (GOO, G12, G15) or ZKron < 19.2 mag 
(G23), using a target catalogue constructed from SDSS DR7 
(G09, G12 and G15) or VST KiDS (G23) imaging. The orig¬ 
inal GAMA concept is described in Driver et al. (2009), the 
tiling algorithm in Robotham et al. (2010), the input cat¬ 
alogue dehnition in Baldry et al.(2010), the optical/near- 
IR imaging pipeline in Hill et al. (2011), the spectroscopic 
pipeline in Hopkins et al. (2013), and the hrst two data 
releases including a complete analysis of the spectroscopic 
campaign and redshift success, in Driver et al. (2011); and 
Liske et al. (2015) respectively. 

One of the scientihc motivations is to assemble a com¬ 
prehensive hux limited sample of ~221,000 galaxies with 
near-complete, robust, fully-sampled spectroscopic coverage 
and robust panchromatic hux measurements from the UV to 
the far-IR and thereafter apply spectral energy distribution 
analysis codes to derive fundamental quantities (e.g., stellar 
mass, dust mass, opacity, dust temperature, star-format ion 
rates etc). 

In this paper we describe the processing and bulk anal¬ 
ysis of the panchromatic data and our discussion is di¬ 
vided into three key sections. Section 2 outlines the gen¬ 
esis and unique pre-processing of each imaging dataset 
into a common astrometric mosaic for each region in each 
band (referred to hereafter as the GAMA SWarps), i.e., ho¬ 
mogenisation of the data. Section 3 outlines our initial ef¬ 
forts towards combining the various hux measurements from 
EUV to far-IR which include a combination of aperture- 
(and seeing)-matched photometry (SDSS/VIKING), ta¬ 
ble matching (GALEX, SDSS/VIKING, WISE), curve-of- 
growth with automated edge detection (GALEX), and op¬ 
tical motivated far-IR source detection (SDSS, SPIRE, 
PAGS). In Section 4 we demonstrate and test the robust¬ 
ness of the PDR. Einally in Section 5 we provide an empirical 
measurement of the EUV-far-IR (0.1 — 500/im) energy out¬ 
put of the Universe in three volume limited slices centred at 
0.5, 1.5, and 2.5 Gyr in look-back time. Note that by energy 
output we refer to the energy being generated per Mpc^ as 
opposed to the energy howing through a Mpc^ (e.g., Driver 
et al. 2008, 2012; Hill et al. 2010). This is important as the 
former refers to the instantaneous energy production rate of 
the Universe (i.e., the luminosity density), whereas the latter 
is the integrated energy production over all time, including 
the relic GMB photons (e.g., Dominquez et al. 2011). 

Throughout this paper we use i7o=70/i7okm Mpc“^ 
and adopt Um = 0.27 and Ua = 0.73 (Komatsu et al. 2011). 
All magnitudes are reported in the AB system. 


2 PANCHROMATIC DATA GENESIS 

Eig. shows the wavelength grasp of the 21 broad-band 
filters. The response curves represent the combined system 
throughputs, normalised to a peak throughput of 1. Also 
shown as a line (in light grey) is the nearby energy output 


© 2010 RAS, MNRAS QQO,|2pOl 


4 


Driver et al. 


GALEXSDSS VL:-IA WISE 

UKIKT 


PACS SPIRE 


CX 

-C 

J 

o 

tn 


0 

0 

0 

0 



Wavelength (/Am) 


o 

CL 

S 

JC, 


p\ 

u 

4j 

d 

pj 


Figure 1. The 21 broad-band combined system throughput curves colour-coded by facility as indicated. Also shown (light grey line) is 
the recently measured(optical)/predicted(mid and far-IR) cosmic spectral energy distribution derived by Driver et al. (2012). This CSED 
can be thought of as an energy weighted “canonical” galaxy spectral energy distribution and highlights how the GAMA PDR filter set 
samples the key energy regime for nearby and low redshift galaxies. Note filters are scaled to a peak throughput of 1 except UKIRT 
which are scaled to 0.5 for clarity. 


Table 1. The GAMA panchromatic imaging regions 


GAMA region 

SWarp RA centre 

SWarp Dec centre 

SWarp ARA 

SWarp A(5 

G09 

09:00:30 

+00:15:00.0 

19dl5m24s 

7d30ml8s 

G12 

11:59:30 

-00:15:00.0 

19dl5m24s 

7d30ml8s 

G15 

14:29:30 

+00:15:00.0 

19dl5m24s 

7d30ml8s 

G23 

23:00:00 

-32:30:00.0 

14d00m00s 

6d00m00s 


Note: G02 is not included here but will be described in a dedicated release paper. 


from the combined z < 0.1 galaxy population derived from 
optical/near-IR analysis of the GAMA dataset (see Driver 
et al. 2012). This highlights how the various bands are sam¬ 
pling the stellar, polycyclic aromatic hydrocarbons (PAHs), 
warm (temperature ~ 50K) and cool (temperature ~ 20K) 
dust emissions of the low redshift galaxy population (the 
curve is shown for the energy output at z = 0). In this 
section we start the process of constructing individual spec¬ 
tral energy distributions (SEDs) for every object within the 
GAMA main survey. 

The hrst step is to place the diverse data onto a common 
astrometric grid. Table. dehnes the extent of the GAMA 
PDR regions. We then use the Terapix SWarp package (see 
Bert in 2010) to build single image mosaics for each wave¬ 
band and each region (see Hill et al. 2011). The SWarp 
package uses the tangent plane (TAN) World Goordinate 
System (WGS) to create a gnomic tangent-plane projection 
centred on the coordinates shown in Table 1. One might 
argue about the merit of constructing such large SWarped 
images (~ 110 deg^ each or up to 80GB for SDSS/VIKING 
data), however it was decided that this was preferable to 
managing the ~ 1 million non-aligned boundaries across the 
PDR. Taking each facility in turn we now describe the pre¬ 
processing necessary to construct our GAMA SWarps. Note 
that in addition to the native-resolution SWarps (see Ta¬ 
ble]^ we also construct a set of SWarps at a common 3.39^^ 
resolution (i.e., 10 times the VISTA pixel scale) for later use 
in deriving coverage hags and background noise estimations. 


2.1 GALEX MIS, GO and Archive data 

The GALaxy Evolution eXplorer (GALEX, Martin et 
al. 2005) was a medium-class explorer mission operated by 
NASA and launched on April 28^^ 2003. The satellite con¬ 
ducted a number of major surveys and observer motivated 
programs, most notably the all-sky imaging survey (AIS; 
typically 200s integrations per tile) and the medium imaging 
survey (MIS; typically 1500s per tile). The GALEX satel¬ 
lite is built around a 0.5-m telescope with a held-of-view 
of 1.13 deg^, a pixel resolution of 1.5'', and a point-spread 
function FWHM of 4.2" and 5.3" in the FUV (153nm) and 
NUV (230nm) bands respectively (Morrissey et al. 2007). 
Imaging data sampled at 1.5 arcsec from V7 of the GALEX 
pipeline forms the basis for constructing the SWarped im¬ 
ages. At the time of commencement of the GAMA survey 
the GAMA regions contained patchy coverage with GALEX. 
A dedicated programme led by one of us (RJT), was pur¬ 
sued providing further GALEX observations to MIS depth 
(1500s) and completed in April 2013 (using funds raised 
from the GAMA and Herschel-ATLAS Gonsortium to reac¬ 
tivate and extend the GALEX mission). The hnal collated 
data provides near-complete NUV and FUV coverage of the 
four primary GAMA regions. Due to the failure of the FUV 
channel mid-mission, the coverage at FUV in G23 is poor. 
However in G09, G12, and G15, coverage is at the 90% level 
in both bands (of which almost ah is at MIS depth in the 
NUV, and 60 percent is at MIS depth in the FUV; see Sec¬ 
tion 2.6). 
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Figure 2. The distribution of exposure times contributing to the 
final FUV and NUV SWarps. The dotted vertical denotes the 
cutoff below which frames are not used in the final SWarps. 

The analysis of the various GALEX datasets are de¬ 
scribed in detail in Andrae (2014) and summarised in Liske 
et al. (2015), and result in background subtracted intensity 
maps scaled to the common GALEX zero-points (Table |^. 
As the data originate from a variety of sources the expo¬ 
sure time is variable (see Eig. [^. To create our SWarps we 
take all available GALEX data frames with exposure times 
greater than 800s. Within the PDR only GALEX has such 
variable integration times. 

In building the SWarps, a common circular mask (of 
radius 35^) was used to trim the outer ~ 5 per cent of the 
image edges where the data quality degrades due to the vi¬ 
gnetting of the telescope aperture (see Morrissey et al. 2007 
and also Drinkwater et al. 2010 who adopted a similar ra¬ 
dius for the WiggleZ survey). In total we have 150, 137, 175, 
and 22 GALEX pointings in EUV and 167, 175, 176 and 133 
in the NUV for GOO, G12, G15 and G23 respectively. These 
are combined to produce single image SWarps at the EUV 
and NUV native resolution for each region. 

Note that a particular subtlety in building the EUV 
and NUV SWarps is the nature of the sky backgrounds. 
In the EUV the majority of pixels have zero flux (i.e., sky 
values of < 1 photon) and hence the distribution of sky 
pixel-values is highly asymmetrical (i.e., Poissonian). Great 
care was taken by the GAMA GALEX team (MS, RJT, 
EA) to model and remove the backgrounds for each individ¬ 
ual frame appropriately and provide to GAMA background 
subtracted EUV data (see Liske et al 2015, section 4.2 for 
further details). Hence when constructing the EUV SWarps 
the background subtraction option was switched off. Eur- 
thermore care should be taken in further background anal¬ 
ysis of GALEX EUV data by only using mean statistics 
and not median statistics because of the highly asymmetri¬ 
cal background distribution, and ensuring sufficient counts 
within any aperture to derive a robust mean. The NUV data 
has a significant sky signal and is therefore processed with 
the SWarp background subtraction on using a 128 x 128 
pixel mesh (i.e., 198"x 198"). 

2.2 SDSS DR7 

The Sloan Digital Sky Survey (York et al. 2000) provides 
uniform optical imaging of the GOO, G12 and G15 regions 
in ugriz bands at 0.4" pixel resolution with a typical PSE 
EWHM of 1.4" (see Hill et al. 2011, figure 3). As the GAMA 


spectroscopic survey was predicated on the SDSS imaging 
(Baldry et al. 2010) there is by design uniform ugriz cover¬ 
age of the three equatorial GAMA regions (GOO, G12 and 
G15). In due course these regions, along with G23, are being 
surveyed by the KiDS team which will provide both deeper 
(2mag) and higher (x2) spatial resolution data (see de Jong 
et al. 2013). 

Here we re-utilise the large mosaic GAMA SWarps built 
from the Sloan Digital Sky Survey Data Release 7 (Abaza- 
jian et al. 2009) by Hill et al. (2011, see update in Liske 
et al. 2015). In brief this involved the construction of both 
native seeing SWarps and SWarps built from data frames 
convolved to a uniform 2" EWHM. The starting point is 
to download all contributing SDSS frames from the DR7 
database, measure the PSE using PSEex (Bertin 2011), 
renormalise the data to a common zero-point, and produce 
both native seeing and convolved data frames (using FGAUSS 
within HEASoft to produce a common PSE EWHM of 2"). 
We then build SWarps at both the native and convolved res¬ 
olutions from the distinct renormalised data frames. During 
the SWarping process (see Bertin 2010; Hill et al. 2011) the 
sky background is subtracted using a coarse 512 x 512 pixel 
median filter to create a grid which in turn is median filtered 
3x3 before being fitted by a bi-cubic spline to represent the 
background structure. The use of a large initial median fil¬ 
ter is to ensure minimal degradation of the photometry and 
shapes of extended systems. 

G23 lies too far south to be observed by SDSS but along 
with GOO, G12 and G15 are being observed to a uniform 
depth within the KiDS survey. The analysis of the KiDS 
data for GAMA and the preparation of the input catalogue 
for G23 will be presented in Moffett et al. (2015). At the 
present time optical SWarps for G23 do not exist. 


2.3 VISTA VIKING 

The Visible and Infrared Telescope for Astronomy (VISTA, 
Sutherland et al. 2015) is a 4.1m short focal length infrared 
optimised survey telescope located 1.5km from the VLT tele¬ 
scopes at Paranal Observatory. VISTA is owned and oper¬ 
ated by ESO and commenced operations on 11^^ December 
2009. VISTA then entered a five year period of survey opera¬ 
tion to conduct a number of ESO Public Surveys (Arnaboldi 
et al. 2007). One of these surveys, the Vista Kilo-degree IN- 
frared Galaxy Survey (VIKING), will cover 1500 deg^ in two 
contiguous regions located in the north and south Galac¬ 
tic caps plus the GOO region. During the first two years of 
operations the VIKING survey prioritised the GAMA and 
Herschel-ATLAS survey regions. The VIKING survey foot¬ 
print therefore covers all four primary GAMA regions (by 
design), in five pass bands {ZYJHKg) at sub-arsecond reso¬ 
lution to projected 5cr point-source sensitivities of 23.1, 22.3, 
22.1, 21.5, 21.2 AB mag (respectively). 

The near-IR camera (VIRGAM, Dalton et al. 2006) 
consists of 16 Raytheon VIRGO HgGdTe arrays (detectors) 
sampling an instantaneous field-of-view of 0.6 deg^ within 
the 1.65 deg diameter field. In routine operation a set of 
micro-dithered and stacked frames are formed, which are 
referred to as PAW-PRINTS. The on-camera dither sequence 
does not cover the gaps between the detectors and hence a 
sequence of six interleaved PAW-PRINTS is required to pro- 
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duce a contiguous coverage rectangular tile of 1.475 deg x 
1.017 deg. 

Paw-print data from the VISTA telescope is pipeline 
processed (Lewis, Irwin V Bunclark 2010) by the Cambridge 
Astronomy Survey Unit (CASU) to produce astrometrically 
and photometrically calibrated data. This process includes 
flat-helding, bias subtraction, and linearity corrections. The 
PAW-PRINTS are then transmitted to the Wide Field Astron¬ 
omy Unit (WFAU) at the Royal Observatory Edinburgh. 
The WFAU combines the PAW-PRINTS into the tiles which 
are then served to the community through both the ESO 
archive and the UK VISTA Science Archive (VSA). As the 
stacked tile data does not include sky-subtraction, sharp 
discontinuities can be introduced into the tiles. An addi¬ 
tional concern is that the tiles may be constructed from 
PAW-PRINTS taken during signihcantly different seeing con¬ 
ditions. As we wish to both sky-subtract and homogenise the 
point-spread function to allow for aperture-matched pho¬ 
tometry (see Hill et al. 2011), we requested access to all 
the VIKING paw-print data provided to the WEAU from 
CASU, which lay within the GAMA primary regions. This 
consisted of 9269 Rice compressed multi-extension fits hies 
(vl.3 data from the CASU archive). These data were ex¬ 
panded out as individual detectors resulting in 148304 indi¬ 
vidual frames. Properties were extracted from the headers 
for each detector (airmass, extinction, exposure time, zero 
point, sky level, seeing) and the seeing measured directly 
using PSEex (Bertin 2011). The data for each individual de¬ 
tector were then rescaled to a common zero point (30) using 
Eqn. 1: 

_ jQ^^]^Q(-0-4(X-2.51ogio(l/t)-(r(secx-l))+Nv.AB-30)) 


where Z is the quoted zero-point, t is the exposure time in 
seconds, r is the extinction in the relevant band and secx 
is the airmass. These values are obtained directly from the 
FITS headers post-CASU processing. Vv.ab is the conversion 
from Vega to AB magnitudes (i.e., 0.521, 0.618, 0.937, 1.384 
or 1.839 for Z,Y,J,H,K respectively) and were derived by 
CASU from the convolution of the complete system response 
functions convolved with the spectrum of Vega and a hat AB 
spectrum. The response functions in comparison to those for 
UKIRT are shown in Eig. 

These data were convolved with the Gaussian kernel re¬ 
quired to produce a EWHM of 2" by assuming the PSE can 
be described as a Gaussian and that the convolution of two 
Gaussians produces a broader Gaussian, i.e., in line with our 
convolved SDSS data (see Hill et al. 2011). Eig.[^shows the 
pre- and post- convolved EWHM as measured by PSEex. 
As can be seen the original seeing is predominantly sub- 
arcsecond as expected from the ESO Paranal (NTT peak) 
site and all the data lies well below our desired target PSE 
EWHM of 2". Because the data is so much better than the 
target PSE EWHM value the assumption of a Gaussian pro- 
hle should produce near-Gaussian hnal PSEs. Note that the 
J band data is observed twice, increasing the abundance 
of independent measurements. The post-processed data is 
centred close to the target PSE EWHM of 2" with some 
indication of slight systematics between the bands at the 5 
per cent level. Note this is not a major concern as we use 
apertures with minimum major or minor diameters of b" 
when measuring our u — Kg aperture-matched photometry 
(see Section 3.1). 



0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 


Seeing (FWHM, arcseconds) 

Figure 3. Pre and post convolution seeing measurements of the 
148304 VISTA VIKING frames using PSFex. 


Erom our initial SWarps we noted that a portion of 
data is clearly of very low quality (see Eig. [^. We there¬ 
fore elected to inspect a subset of the data by selecting 
three categories: OUTLIERS dehned as those with seeing bet¬ 
ter than 0.5^^ or worse than 1.5^\ a zero-point multiplier of 
greater than 40, a sky value of less than 100 ADU counts or 
a GASU TILECODE not equal to 0, 56, or -1, i.e., 9535 frames 
in total; CONTROL defined as a random set of 1000 frames 
not included in the above selection; and SPARSE defined as 
every detector 8 frame not already included in one of the 
earlier samples, i.e., 6945 frames. These 16590 frames were 
inspected by two of us (SPD, AHW) using the MOGRIFY rou¬ 
tine within the IMAGEMAGICK package to generate greyscale 
images where the lowest 2 per cent of data were set black, 
the highest 10 per cent white, and with histogram equalisa¬ 
tion in-between. This scaling amplifies background gradients 
rendering even the best quality data in the poorest light (see 
Fig.Q. We then rejected or accepted the frames via visual 
inspection and attempted to identify a measurable quantity 
which best separated out the rejected frames, see Eig. 
This resulted in the adoption of a simple cut on the zero- 
point multiplier factor, whereby all frames which require a 
rescaling of x30 or more are rejected in addition to those 
already identihed from the visual inspections. In total 3262 
of our 148304 frames were rejected (i.e., 2.2 per cent of the 
data). Examples of accepted and rejected frames are shown 
in Eig. and common causes are bright sky, detector failures 
and telescope pointing errors. 

The remaining frames were then SWarped (Bertin et 
al. 2010) to the GAMA PDR regions specified in Table 
with a pixel size of 0.339" using the TAN WGS projection. 
During the SWarping process the background for each con¬ 
tributing detector was removed using a 128 x 128 pixel me¬ 
dian hlter which in turn was median hltered by a 3 x 3 
grid before being fitted with a bi-cubic spline. The choice of 
background hlter size is critical; too high and the structure 
of the tiling becomes apparent in the SWarp (see Eig. 6), 
too low and galaxy photometry can be affected (see Eig. 7). 
To optimise the background hlter size we produced frames 
with a range of background hlter sizes and performed struc¬ 
tural analysis of the brightest 100 galaxies using SIGMA 
(Kelvin et al. 2012). Eig.[^shows the magnitude offsets and 
Eig. E shows how the measured major-axis half-light radii 
vary with background mesh size. We tested pixel grids of 
512 X 512, 256 x 256, 128 x 128 and 64 x 64 and only the 
smallest hlter size had any noticeable impact on the mea¬ 
sured properties and hence the second smallest hlter size was 
adopted. Note that this hner hltering (compared to SDSS) 
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Figure 4. Examples of poor quality (top three rows) and accept¬ 
able quality (bottom row) VIKING frames. Approximately 12 per 
cent of the VIKING data were visually inspected based on out¬ 
lying values in airmass, sky background, zero-point, and seeing. 
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Zeropoint multiplier 


100 


Figure 5. Seeing versus zero-point multiplier for each band (as 
indicated). A cut of 30 appears to isolate the majority of low 
quality frames (indicated by the coloured points). 



Figure 7. An illustration of the impact of over-smoothing the 
background on extended objects. The galaxy shown is the largest 
system in the GAMA region, NGG0895 (located in G02). In the 
rightmost panel a significant portion of the galaxy has been re¬ 
moved due to the 64 x 64 pixel sky-subtraction process. 


is absolutely necessary because of (a) the mode of observa¬ 
tion (pointed v drift-scan) and (b) the higher-degree of sky 
spatial variations in the near-IR wavebands. 

2.4 WISE 

The Wide-Field Infrared Survey Explorer (WISE; Wright 
et al. 2010) is a medium-class explorer mission operated by 
NASA and was launched on 14^^ December 2009. Eollowing 
approximately 1 month of checks WISE completed a shallow 
survey of the entire sky in 4 infrared bands (3.4, 4.6, 12 and 
22/im) over a ten month period. WISE is built around a 40- 
cm telescope with a 47' x 47' field-of-view, and scans the sky 
with an effective exposure time of 11s per frame. Each region 
of sky is typically scanned from tens to hundreds of times 
(with fields further from the ecliptic being observed more 
frequently). This allows the construction of deep stacked 
frames reaching a minimum 5cr point source sensitivity of 
0.08, 0.11, 0.8 and 4 mJy in the Wl(3.4/im), W2(4.6//m), 
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Magnitude (1 Ore) 



Magnitude (lOre) 



Magnitude (lOre) 


Figure 8. A comparison of flux measurements of the brightest 
100 galaxies with varying background subtraction meshes. In each 
case the flux is compared against that derived from the 512 x 512 
pixel background mesh. In general it is only galaxies brighter 
than 14^^ magnitude with the 64 x 64 background mesh whose 
photometry is compromised. 


W3(12/im) and W4(22/im) bands (see Wright et al. 2010). 
The base “Atlas” data consists of direct stacks and associ¬ 
ated source catalogues which are publicly available via the 
WISE and AllWISE data release hosted by the Infrared Sci¬ 
ence Archive (IRSA). These public data have point-spread 
function EWHM resolutions of ~ 8.4", 9.2", 11.4" and 18.6" 
in Wl, W2, W3 and W4 respectively and a 1.375"/pixel 
scale. However, because of the stability of the point-spread 
function of the WISE system, higher resolution can be at¬ 
tained using deconvolution techniques, in particular “driz¬ 
zled” co-addition and the Maximum Correlation Method of 
Masci & Eowler (2009; see Jarrett et al. 2012). Here we use 
data which has been re-stacked via the drizzle method as 
the MCM or HiRes method is computationally expensive 
and only suited for very large nearby galaxies (see Jarrett 
et al. 2012, 2013). In brief this involves: 

(1) gain-matching and rescaling the data ensuring a common 
photometric zero-point calibration, 

(2) background level offset-matching, 

(3) flagging and outlier rejection. 


Half-Light Radius 



Sersic Index 



Sersic index 


Figure 9. The major-axis half-light radius (upper panel) and 
the Sersic index distributions (lower panel) for the brightest 100 
galaxies for various background mesh sizes as indicated. 


(4) co-addition using overlap area weighted interpolation 
and drizzle. 

Here drizzle refers to the Variable Pixel Linear Re¬ 
construction technique of co-addition using a Point Re¬ 
sponse Function kernel to construct the mosaics. Full details 
are provided in the WISE ICORE documentation (Masci 
2013). The drizzled data results in final point-source EWHM 
of 5.9", 6.5", 7.0", and 12.4" (respectively), see Cluver et 


shows a comparison for one of our GAMA galaxies between 
the “Atlas” and drizzled image in each of the four bands. The 
“drizzled” frames are provided to the GAMA team stacked, 
calibrated to a common zero point, and background sub¬ 
tracted in sections of 1.56° x 1.56°. These frames are then 
SWarped into a single large mosaic at the native pixel res¬ 
olution of 1" using the same field centre, and projection 
system as for the previous datasets. In re-gridding the data 
we also include the SWARP background subtraction using 
a 256 X 256 pixel filter. 


al. (2014) and Jarrett et al. (2012) for further details. Fig. 10 


2.5 Herschel-ATLAS 

The Herschel Space Observatory (Pilbratt et al. 2010) is 
operated by the European Space Agency and was launched 
on May 14^^ 2009 and conducted a number of major sur¬ 
vey campaigns during its 3.5 years of operation. The largest 
extragalactic survey, in terms of areal coverage, is The Her¬ 
schel Astrophysical Terahertz Large Area Survey (Herschel- 
ATLAS; Eales et al. 2010). Herschel-ATLAS images were ob¬ 
tained using Herschel’s fast-scan parallel mode and covered 
~600 deg^ of sky in five distinct sky regions which included 
the four principal GAMA fields. The co-ordinated observa¬ 
tions used both the PAGS (Poglitsch et al. 2010) and SPIRE 
(Griffin et al. 2010) instruments to obtain scans at 100/im, 
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Figure 10. A visual comparison of WISE Atlas images (upper rows) and WISE drizzled images (lower rows) for bands WI, W2, W3, 
W4 (left-to-right). The panels are displayed over comparable ranges and the improvement in spatial resolution via the drizzling technique 
is self-evident. For details of the drizzling technique see Jarrett et al. (2012). 


160/im, 250/im, 350/im, 500/im, i.e., sampling the warm and 
cold dust components of galaxies from z=0 to z=4. The fi¬ 
nal maps were the combination of two orthogonal cross-scans 
giving rise to PSFs with Gaussianised FWHM of 9.6'^ and 
12.5" in 100 and IGOfim and 18", 25" and 35" in the 250, 
350 and 500 iim bands respectively (see Valiante et al. 2015 
for full details of the PSF characterisation). The data were 
processed, calibrated, nebularized to remove large scale fluc¬ 
tuations due to cirrus and large scale clustering of high-z 
sources, (see Valiante et al.2015. and Maddox et al. 2015), 
and finally mosaiced by the Herschel-ATLAS data reduction 
team who provided the final maps and 5cr source detection 
catalogues. The reduction process for the two instruments 
are described in detail in Ibar et al. (2010), and Pascale et 
al. (2011), to be superseded shortly by Valiante et al. (2015), 
and the method for source detection is described in detail 
in Rigby et al. (2011) also updated in Valiante et al. (2015). 
The absolute zero point calibration is accurate to ±10 per 
cent for PACS and ±7 per cent for SPIRE which provides a 
potential systematic pedestal in addition to the random sky 
and object photon noise errors estimated later. Note that 
as the Herschel-ATLAS data have not yet been publicly re¬ 
leased they remain subject to change. Every attempt will 
be made to ensure the online GAMA PDR provides notifi¬ 
cations of any changes or updates. 

To date the Herschel-ATLAS data have been used to 
study the dust and star-formation properties of both near 
and distant galaxies based on far-IR/optical matched sam¬ 
ples (see for example Dunne et al. 2011; Smith et al. 2011, 
2012; Bourne et al. 2012 Rowlands et al. 2012). To pre¬ 
prepare the data for GAMA we re-SWarp the mosaics pro¬ 
vided onto a uniform grid using the field centres from Ta¬ 
ble using the TAN WGS projection, and preserving the 
original pixel size as specified in the file headers and shown 
on Tabled 


2.6 Cosmetic and noise characteristics of the 
GAMA SWarp set 

To assess the quality of GAMA SWarps we derive the 
background noise distributions (i.e., sky-subtracted), within 
selected regions for each of our LOW-RES (i.e., 3.39") 
SWarps, which are displayed from —2cr to ±2cr in Eigs. 43 
to A|^for G09, G12, G15 and G23 respectively. The black 
rectangle represents the GAMA region and the dotted blue 
rectangle the selected region from which the noise charac¬ 
teristics are derived (the mode and 3cr-clipped standard de¬ 
viation). These images show no obvious major sky gradients 
across the sky regions, however, they do show interesting 
substructure which highlights correlations in the underly¬ 
ing noise properties. In most cases the correlations high¬ 
light the genesis, i.e., the SDSS stripes, GALEX pointings, 
and VIKING paw-prints. In these cases the noise properties 
for each particular frame/scan is dictated by the conditions 
during observations (SDSS and VISTA) or the variability 
of the various integration times (GALEX). While uniform 
backgrounds are highly desirable, these are never achieved 
in practice. Some SDSS scans will be slightly less noisy than 
others and some paw-prints will have signihcantly amplihed 
noise characteristics. Interestingly the WISE and Herschel- 
ATLAS data show the least structure which mainly rehects 
the benehts of using hxed integration times as well as op¬ 
erating outside the conhnes of a time-varying atmosphere. 
However, some impact of observing close to the moon is 
apparent in the WISE G12 SWarps. Also noticeable in the 
Herschel-ATLAS data is the reduced noise in the overlap 
regions as expected. 

The noise distributions derived from the GAMA 
SWarps are shown in Table |3] for GALEX, SDSS, VISTA 
and WISE data these statistics are derived from htting a 
Gaussian distribution to the histogram of data values below 
the mode. They therefore do not include any confusion esti- 
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Table 2. Key meta-data information of the contributing datasets. 


Facility 

Dataset 

Instrument 

Filter 

Pivot 

Pixel 

Point-source 

Frames 

”*AB — ”*Vega 


or survey 

or technique 

name 

Wavelength 

Resolution 

FWHM 

Supplied 

(mag) 

GALEX 

MIS+GO 

- 

FUV 

1535A 

1.5" 

4.1" 

279 

2.16 

GALEX 

MIS+GO 

- 

NUV 

2301A 

1.5" 

5.2" 

297 

1.67 

SDSS 

DR7 

- 

u 

3557A 

0.339" 

1.4" 

26758 

0.98 

SDSS 

DR7 

- 

9 

4702A 

0.339" 

1.4" 

26758 

-0.10 

SDSS 

DR7 

- 

r 

6175A 

0.339" 

1.4" 

26758 

0.15 

SDSS 

DR7 

- 

i 

7491A 

0.339" 

1.4" 

26758 

0.38 

SDSS 

DR7 

- 

z 

8946A 

0.339" 

1.4" 

26758 

0.54 

VISTA 

VIKING 

VIRGAM 

Z 

88005 

0.339" 

0.85" 

15360 

0.521 

VISTA 

VIKING 

VIRGAM 

Y 

10213A 

0.339" 

0.85" 

15797 

0.618 

VISTA 

VIKING 

VIRGAM 

J 

12525A 

0.339" 

0.85" 

34076 

0.937 

VISTA 

VIKING 

VIRGAM 

H 

16433A 

0.339" 

0.85" 

15551 

1.384 

VISTA 

VIKING 

VIRGAM 

Ks 

21503A 

0.339" 

0.85" 

16340 

1.839 

WISE 

AllSky 

drizzled 

W1 

3.37/rm 

1 " 

5.9" 

40 

2.683 

WISE 

AllSky 

drizzled 

W2 

4.62/rm 

1" 

6.5" 

40 

3.319 

WISE 

AllSky 

drizzled 

VF3 



7.0" 

40 

5.242 

WISE 

AllSky 

drizzled 

VF4 

22.8/rm 

1 " 

12.4" 

40 

7.871 

Herschel 

ATLAS 

PAGS 

lOOfim 

101p,m 

3" 

9.6" 

4 (& 1 for G23) 

N/A 

Herschel 

ATLAS 

PAGS 

1Q0/J,m 

161/rm 

A" 

12.5" 

4 (& 1 for G23) 

N/A 

Herschel 

ATLAS 

SPIRE 

250/rm 

249p,m 

6" 

18" 

4 (& 1 for G23) 

N/A 

Herschel 

ATLAS 

SPIRE 

350/rm 

357/rm 

8" 

25" 

4 (& 1 for G23) 

N/A 

Herschel 

ATLAS 

SPIRE 

500p,m 

504p,m 

12 " 

36" 

4 (& 1 for G23) 

N/A 
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Figure 11. The sensitivity reached in each band as derived from the GAMA SWarps (black lines) and compared to the listed values 
(coloured lines). Also shown in grey is a typical SED for a dusty galaxy with tab = 19-8 mag. 


mate and assume the noise is uncorrelated. In all cases the 
distributions are very well described by a normal distribu¬ 
tion implying that the systematic frame-pistoning (i.e., ZP 
offsets) in the data (arising from the independent calibra¬ 
tion of the distinct pointings), is operating at a relatively 
low level and within the range of the pixel-to-pixel varia¬ 
tions. Using the 3cr-clipped standard deviations we derive 
(analytically) the Icr surface brightness limits and the 5cr 
point-source detection limits for each of the SWarp images 
(see Table. [^. For the PACS and SPIRE data, where cor¬ 
related noise is believed to be an issue, we derive the 5cr 
detection limits directly by placing apertures equivalent to 


the Beam size at random locations across the SWarps and 
measuring the standard deviation of the resulting aperture 
fluxes (again fitting to the distribution below the mode). 
In Fig. m the GAMA SWarp detection limits are com¬ 
pared to the values listed online for each facility (as indi¬ 
cated by the colour lines). For GALEX MIS, SDSS DR7, 
and VIKING, the depths probed agree extremely well. Note 
that our derived WISE W1 band limit appears significantly 
deeper than that quoted by the WISE collaboration this is 
because our values ignore confusion (i.e., fits to the negative 
noise distribution) whereas the WISE quoted value incorpo¬ 
rates this aspect. For Herschel-ATLAS SPIRE we note that 
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Table 3. Surface brightness limits of our GAMA SWarp set (FUV to mid-IR). 


SWarp 

Facility/Filter/Field 

Zero-Point 
(AB mag for lADU) 

SWarp mean 
(/") 

Ic^Sky 

(mag arcsec“^)'t' 

5(7 limit 

(mag)t (Jy)^ 

Goverage 

(%) 

GALEX FUV G09 

18.82 

0.000148 

28.41 

25.23 

2.94E - 07 

88 

GALEX FUV G12 

18.82 

5.04E - 05 

29.58 

26.40 

IE-07 

92 

GALEX FUV G15 

18.82 

0.00013 

28.54 

25.37 

2.59E - 07 

95 

GALEX FUV G23 

18.82 

0.000266 

27.77 

24.59 

5.31E-07 

75 

GALEX NUV G09 

20.08 

0.00125 

27.35 

23.92 

9.84E - 07 

94 

GALEX NUV G12 

20.08 

0.00116 

27.43 

23.99 

9.17E-07 

97 

GALEX NUV G15 

20.08 

0.00161 

27.07 

23.64 

1.27E - 06 

95 

GALEX NUV G23 

20.08 

0.00109 

27.50 

24.07 

8.58E - 07 

99 

SDSS u G09 

30.00 

142 

24.61 

22.24 

4.61E - 06 

100 

SDSS u G12 

30.00 

163 

24.47 

22.09 

5.28E - 06 

100 

SDSS u G15 

30.00 

156 

24.51 

22.14 

5.06E - 06 

100 

SDSS g G09 

30.00 

54.3 

25.66 

23.29 

1.76E - 06 

100 

SDSS g G12 

30.00 

65.3 

25.46 

23.09 

2.12E - 06 

100 

SDSS g G15 

30.00 

62.6 

25.50 

23.13 

2.03E - 06 

100 

SDSS r G09 

30.00 

76.7 

25.28 

22.91 

2.49E - 06 

100 

SDSS r G12 

30.00 

95.4 

25.05 

22.67 

3.09E - 06 

100 

SDSS r G15 

30.00 

90.5 

25.10 

22.73 

2.93E - 06 

100 

SDSS i G09 

30.00 

116 

24.84 

22.47 

3.75E - 06 

100 

SDSS i G12 

30.00 

140 

24.62 

22.25 

4.56E - 06 

100 

SDSS i G15 

30.00 

129 

24.72 

22.35 

4.19E-06 

100 

SDSS z G09 

30.00 

506 

23.23 

20.86 

1.65E - 05 

100 

SDSS ^ G12 

30.00 

579 

23.09 

20.71 

1.88E - 05 

100 

SDSS z G15 

30.00 

556 

23.13 

20.76 

1.81E - 05 

100 

VIKING Z G09 

30.00 

59.8 

25.55 

23.18 

1.94E - 06 

100 

VIKING Z G12 

30.00 

60.6 

25.54 

23.17 

1.97E - 06 

100 

VIKING Z G15 

30.00 

62 

25.51 

23.14 

2.01E - 06 

99 

VIKING Z G23 

30.00 

67.9 

25.41 

23.04 

2.2E - 06 

100 

VIKING Y G09 

30.00 

123 

24.77 

22.40 

3.98E - 06 

100 

VIKING Y G12 

30.00 

110 

24.89 

22.52 

3.56E - 06 

100 

VIKING Y G15 

30.00 

111 

24.88 

22.51 

3.61E-06 

100 

VIKING Y G23 

30.00 

129 

24.71 

22.34 

4.2E - 06 

100 

VIKING J G09 

30.00 

167 

24.44 

22.06 

5.43E - 06 

100 

VIKING J G12 

30.00 

161 

24.48 

22.10 

5.23E - 06 

100 

VIKING J G15 

30.00 

146 

24.58 

22.21 

4.74E - 06 

100 

VIKING J G23 

30.00 

155 

24.52 

22.14 

5.04E - 06 

100 

VIKING H G09 

30.00 

329 

23.70 

21.33 

1.07E - 05 

98 

VIKING H G12 

30.00 

302 

23.79 

21.42 

9.82E - 06 

99 

VIKING H G15 

30.00 

313 

23.75 

21.38 

1.02E - 05 

97 

VIKING H G23 

30.00 

325 

23.71 

21.34 

1.06E - 05 

100 

VIKING K G09 

30.00 

332 

23.69 

21.32 

1.08E - 05 

100 

VIKING K G12 

30.00 

337 

23.67 

21.30 

1.09E - 05 

100 

VIKING K G15 

30.00 

303 

23.79 

21.42 

9.83E - 06 

100 

VIKING K G23 

30.00 

285 

23.86 

21.48 

9.25E - 06 

100 

WISE W1 G09 

23.18 

0.262 

24.64 

21.09 

1.33E - 05 

100 

WISE W1 G12 

23.18 

0.281 

24.56 

21.01 

1.43E - 05 

100 

WISE W1 G15 

23.14 

0.21 

24.84 

21.29 

l.llE-05 

100 

WISE W1 G23 

23.14 

0.187 

24.96 

21.41 

9.9E - 06 

100 

WISE W2 G09 

22.82 

0.327 

24.04 

20.38 

2.55E - 05 

100 

WISE W2 G12 

22.82 

0.367 

23.91 

20.26 

2.87E - 05 

100 

WISE W2 G15 

22.82 

0.264 

24.27 

20.61 

2.06E - 05 

100 

WISE W2 G23 

22.82 

0.229 

24.42 

20.77 

1.79E - 05 

100 


^ lAia = ZP - 2.5 \ogiQ(aAnij)- 

^ 5(7 limit = ZP — 2.5 log]^Q(5V ttHWHM"^ a adu) where HWHM is Half Width Half-Maximum of the seeing-disc (i.e., 0.5 FWHM). 

o Fjy{Jy) = 3631 X 


the agreement with the SPIRE values reported in Valiante 
et al. (2015) is extremely good. 


2.7 Astrometric verification 

To check the astrometric alignment of the SWarps we 
run SExtractor over either the entire SWarp set (GALEX, 


WISE, and Herschel) or a 4sq deg section from G12 (SDSS 
and VIKING, to keep GPU requirement manageable). We 
use relatively high signal-to-noise cuts of 100 (GALEX, 
where the data is not uniform), 10 (WISE and Herschel), or 
3 (SDSS and VIKING). We then match to either the GAMA 
InputGat with r < 17.0 mag (GALEX, SDSS, VIKING 
and WISE) or the GAMA TilingGatv45 with r < 17.0 mag 
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Table 4. Surface brightness limits of our GAMA SWarp set (mid-IR to far-IR. 


SWarp 

Facility/Filter/Field 

Zero-Point 
(AB mag for lADU) 

SWarp mean 
(/") 

Ic^Sky 

(mag arcsec“^)'t' 

5f7 limit 

(mag)t (Jy)^ 

Goverage 

(%) 

WISE W3 G09 

23.24 

2.33 

22.32 

18.59 

0.000133 

100 

WISE W3 G12 

23.24 

2.68 

22.17 

18.44 

0.000153 

100 

WISE W3 G15 

23.24 

1.77 

22.62 

18.89 

0.000101 

100 

WISE W3 G23 

23.24 

2.18 

22.39 

18.66 

0.000125 

100 

WISE W4 G09 

19.60 

0.278 

20.99 

16.64 

0.000802 

100 

WISE W4 G12 

19.60 

0.305 

20.89 

16.54 

0.000879 

100 

WISE W4 G15 

19.60 

0.208 

21.31 

16.96 

0.000599 

100 

WISE W4 G23 

19.60 

0.265 

21.05 

16.69 

0.000762 

100 

PACS 100 G09 

8.90 

0.000562 

- 

12.96 

0.0894 

100 

PACS 100 G12 

8.90 

0.000545 

- 

12.99 

0.0879 

100 

PACS 100 G15 

8.90 

0.000547 

- 

12.99 

0.0863 

100 

PACS 100 G23 

8.90 

0.000476 

- 

13.14 

0.0795 

100 

PACS 160 G09 

8.90 

0.000278 

- 

13.44 

0.103 

100 

PACS 160 G12 

8.90 

0.000273 

- 

13.46 

0.101 

100 

PACS 160 G15 

8.90 

0.000271 

- 

13.47 

0.101 

100 

PACS 160 G23 

8.90 

0.000227 

- 

13.66 

0.0903 

100 

SPIRE'’ 250 009 

11.68 

0.000759 

- 

12.56 

0.0343 

80 

SPIRE'’ 250 G12 

11.68 

0.00073 

- 

12.60 

0.0330 

81 

SPIRE^ 250 G15 

11.68 

0.00073 

- 

12.59 

0.0333 

84 

spire'’ 250 G23 

11.68 

0.000885 

- 

12.52 

0.0357 

100 

SPIRE^ 350 G09 

11.67 

0.000447 

- 

12.36 

0.0412 

80 

spire'’ 350 G12 

11.67 

0.000424 

- 

12.41 

0.0394 

81 

spire'’ 350 G15 

11.67 

0.000423 

- 

12.41 

0.0393 

84 

SPIRE'’ 350 G23 

11.67 

0.000518 

- 

12.51 

0.0357 

100 

SPIRE'’ 500 009 

11.62 

0.000228 

- 

12.16 

0.0495 

80 

SPIRE^ 500 G12 

11.62 

0.000217 

- 

12.23 

0.0467 

81 

spire'’ 500 G15 

11.62 

0.000221 

- 

12.21 

0.0476 

84 

SPIRE^ 500 G23 

11.62 

0.000257 

- 

12.17 

0.0490 

100 


Micr = - 2.51ogio(o-ADu)- 

^ 5cr limit = ZP — 2.5 log]^Q(5\/ HWHM‘^a where HWHM is Half Width Half-Maximum of the seeing-disc (i.e., 0.5 FWHM). 

o F^(Jy) = 3631 X 

^ SPIRE maps are in units of Jansky per Beam and to generate these zero-points we have added a factor 2.5 log]^Q(H/A^) where B is 
the beam size given as 466, 821, and 1770 sq arcsec in 250, 350 and 500/rm respectively and N is the pixel size given in Table 


and {u — g) < 1.5 (Herschel-ATLAS) to isolate star-forming 
galaxies. Fig, shows the resulting A RA and A Dec dia¬ 
grams for each band compared to the canonical r-band data 
(grey data points). On Fig. 12 the blue cross (mostly not vis¬ 
ible) defines the centroid and the thick green circle indicates 
the PSF FWHM for that band. The thick blue band defines 
the region which encloses 66 percent of the population (after 
accounting for the density of random mis-matches), and the 
thin blue circles enclose either 50 percent or 80 percent of 
valid matches. Fig. highlights that in ah cases the cen¬ 
troid of the RA and Dec offset is extremely close to zero 
(below 0.3'' in ah bands with the FUV and NUV showing 
the largest offsets, and below 0.02" in the optical and near- 
IR), and that the 66 percent sprawl lies within 0.5x the PSF 
FWHM in ah bands. We therefore consider the astrometry 
to be as one would expect given the respective FWHM see¬ 
ing values. 


2.8 Visual inspection of the combined data and 
data access 

Our full dataset is diverse and the volume large. In order 
to inspect the data we have developed a publicly available 
online tool which provides both download links to the indi¬ 
vidual SWarps, as well as an option to extract image regions 
from the dataset. Users can also build RGB colour images 




Figure 13. A comparison of the quality of the SDSS 2 ; band data 
(left) against the VISTA VIKING Z band (right). 


using any of the 21 bands as well as overlay contours and 
basic catalogue information (e.g., GAMA IDs, photometry 
apertures, and object locations). The GAMA Panchromatic 
SWarp Imager (4^) is therefore extremely versatile and use¬ 
ful for exploring the data volume: 
http://gama-psi.icrar.org/ 

Fig. show examples of various extractions us¬ 

ing the tool with Fig. |13| showing the signihcant increase in 
depth from the SDSS z band data to the VISTA VIKING Z 
band. Fig, p^ shows a single GAMA galaxy in 20 of the 21 
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Figure 12. Confirmation of the astrometric accuracy. Each diagram shows the positional offsets of that particular band against the r 
band GAMA Input or Tiling catalogues (grey data points). The centroid of the population is shown with a blue cross and the (native) 
PSF FWHM is shown as a green circle. The circles enclosing 50, 66 and 80 per cent of the population is shown by a thin, thick and thin 
blue line respectively. In all cases the relative astrometry is robust to < O.IPSF FWHM, and the 66 per cent spread enclosed with 0.5x 
the PSF FWHM. 


bands (note that the SDSS z band is not shown here), and 
Fig. shows various colour combinations with contours, 
IDs and apertures overlaid as indicated. Note that searches 
can be made based on GAMA ID or RA and Dec and is 
therefore of use to high-z teams with objects in the GAMA 
regions (e.g., Herschel-ATLAS team). 

Using GAMA ^ via the link above, one can also ac¬ 
cess the individual SWarps files including the native, con¬ 
volved and weight-maps and the XML files which contain, 
the pixel data, a description of the weights, and a listing of 
the constituent files making up the SWarp respectively. The 
weight-maps are particularly useful and can be used to de¬ 
termine both the coverage and provide a mask. Zero values 
in the weight SWarp imply no data while non-zero values 
imply coverage. These weight-maps have been used to gen¬ 
erate the coverage statistics shown in col. 7 of Table. The 
SWarps, weight-maps, and XML files can all be downloaded 
from: 

http://gama-psi.icrar.org/panchromaticDR.php 


However, note that files sizes vary from 100KB (for XML 
files) to up to 80GB (for SWarps and weight-maps). 


3 PANCHROMATIC PHOTOMETRY FOR G09, 
G12 AND G15 ONLY 

Vital to successful analysis of panchromatic data are robust 
flux measurements, robust errors, and a common deblending 
solution. This is particularly difficult when the flux sensitiv¬ 
ities and spatial resolutions vary significantly, as is the case 
with the GAMA PDR (see Figs[^[^&[^ i.e., 35" to 0.7" 
spatial resolution). In an ideal situation one would define 
an aperture in a single band and then place the same aper¬ 
ture at the same astrometric location in data with identical 
spatial sampling. This is the strategy we pursued in Hill 
et al. (2011, see also Driver et al. 2011) to derive u — Kg 
aperture-matched photometry (using the seeing-convolved 
SWarps convolved to a 2" FWHM). While we can still im¬ 
plement this strategy in the u to Ks range (see Section 4.2 
below) we cannot easily extend it outside this wavelength 
range because of the severe resolution mismatch (see Ta- 
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GALEX FUV 




SDSSi 






Herschel S350 




SDSS 


VIKING Y 

f% ■ 

WISE 1 


Herschefpiod 


Herschel S500 


Figure 14. 20 band panchromatic imaging for a 1' x 1' region 
centred on GAMA galaxy G371633. Filters increase in wavelength 
proceeding from left to right and top to bottom, note the SDSS 
^ filter is omitted, i.e., FUV, NUV, u, g, r, i, Z, Y, J, H, K, Wl, 
W2, W3, W4, 100, 160, 250, 350, 500. Produced using the GAMA 
4^: http://gama-psi.icrar.org/ 



217.064 217.051 217.037 217.024 

RA'd9g 


217.064 217.051 217.037 217.024 

RA'dgg 


Figure 15. Upper left: A colour composite image of G79334 pro¬ 
duced by combining the SDSS g Sz, r with the VIKING H band 
images. Overlaid are the contours from SPIRE-250 band. GAMA 
IDs are marked. Upper right: A colour composite of G79334 from 
GALEX NUV, VIKING Z & Ks and with contours overlaid from 
WISE. Lower left: A 2' x 2' colour composite centred on G48432 
made from data extracted from the GAMA SDSS r Sz i SWarp 
combined with the VIKING H SWarp and with the apertures 
used for the aperture-matched photometry overlaid, (lower right) 
A composite colour image of G48432 made from GALEX FUV 
(blue channel), WISE Wl (green channel) and SPIRE 250/j,m 
(red channel and contours). All images produced using the online 
GAMA ^ tool: http://gama-psi.icrar.org/ 


ble[^. Software (LAMBDAR) is being developed to specif¬ 
ically address this issues and will be described in a compan¬ 
ion paper (Wright et al. 2015). In the meantime, we assem¬ 
ble a benchmark panchromatic catalogue from a combina¬ 
tion of aperture-matched photometry, table matching, and 
optically-motivated (forced) photometry. It is worth noting 
that the GAMA PDR assembled here while heterogeneous 
across facilities is essentially optimised for each facility, and 
therefore optimal for studies not requiring broad panchro¬ 
matic coverage. 

In the FUV and NUV, we use the GAMA GALEX cat¬ 
alogue described in Liske et al. (2015) and which uses a va¬ 
riety of photometry measures including curve-of-growth and 
the GALEX pipeline fluxes. In the optical and near-IR we 
apply the aperture-matched method mentioned above and 
described in detail in the next sections. In the mid-IR we use 
the WISE catalogues described in Gluver et al. (2014). Long- 
wards of the WISE bands we adopt a strategy developed by 
the Herschel-ATLAS team (Bourne et al. 2012, see Appendix 
A) to produce optically-motivated aperture measurements 
(sometimes referred to as forced photometry). This is ap¬ 
plied to all GAMA targets which lie within the PAGS and 
SPIRE 100 to 500/im data. 


3.1 Aperture-matched photometry from u to Kgi 
IOTA 

The u to Ks band data has been convolved to a common 2" 
EWHM seeing (see Eig. [^. Eor each object in the GAMA 
tiling catalogue with a secure redshift (TilingGatv44, i.e., a 
valid galaxy target within the specified regions with tab < 
19.8 mag, see Baldry et al. 2010) we perform the following 
tasks: 

(1) extract a 1001 x 1001 pixel region in all 10 bands 
{ugrizZYJHKs), 

(2) run SExtractor in dual object mode with r as the pri¬ 
mary band, 

(3) identify the SExtractor object closest to the central pixel 
(2'' max), 

(4) extract the photometry for this object in the two bands, 

(5) repeat for all bands. 

In essence this process relies on SDSS DR7 for the ini¬ 
tial source detection and initial classification including an 
r-band Petrosian flux limit to define the input catalogue. 
However, the final deblending and photometry is ultimately 
based on SExtractor (using the parameters described in 
Liske et al. 2015 optimised for our convolved data). An iden¬ 
tical aperture and mask and deblend solution — initially de¬ 
fined in the r band — is then applied to the ugizZY JHKg 
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bands. In order to manage this process efficiently for 220k 
objects we use an in-house software wrapper, IOTA. 


3.2 Recalibration of the u to Kg photometry 


The VIKING data is relatively new and to assess the abso¬ 
lute zero-point errors we test the consistency of the photom¬ 
etry between our measured VIKING data and the 2MASS 
point source catalogue. We achieve this by extracting all 
catalogued stars in the extended GAMA regions from Input- 
GatvOG which itself is derived from SDSS DR7 (see Baldry 
et al. 2010). To obtain near-IR flux measurements we up¬ 
loaded the objects classified as stars (see Baldry et al. 2010) 
to the IPAG Infrared Science Archive (IRSA) and queried 
the 2MASS All-Sky Point Source Gatalogue (on 2013-06-07). 
We obtained 498,637 matches for which photometry existed 
in 1 or more of the 2MASS bands {JHKg). This sample was 
trimmed to the exact GAMA RA extents to produce cata¬ 
logues of 201671, 92224 and 131976 stars in G09, G12, and 
G15 respectively. We ran IOTA on these objects to derive 
ugrizZY JHKs photometry based on Kron a pert ures with 
a minimum aperture diameter of b". Fig. 16 shows the 
resulting zero-point comparisons versus magnitude (left pan¬ 
els) and versus the VIKING (J — K)ab colour (right pan¬ 
els) for filters ugrizZ JHKg (top to bottom) respectively. 
Note that for the ugriz bands we compare directly to SDSS 
PSF mags corrected to AB (i.e., uab = usbss — 0.04 and 
ZAB = ^SDSS + 0.02) for the ZJHKg bands we convert the 
2MASS mags into the VISTA passband system, using the 
colour transformations derived by the VISTA Variables in 
the Via Lactea Survey (VVV) team (Soto et al. 2013) which 


TvISTA = e/2MASS — 0.077(J2MASS — 772MASs) (2) 

77vista = 772MASS + 0.032( J2MASS — 772mass) (3) 

i^VISTA = i^2MASS + 0.010( J2MASS — i^2MASs) (4) 

Finally we implement the Vega to AB correction appropriate 
for the VISTA filters, see Table. 

At brighter magnitudes the deeper VIKING data will 
suffer from saturation, and at fainter magnitudes the shal¬ 
low 2MASS data will become swamped by noise. Fig, and 
17 show the direct comparisons (black data points), and the 
data which we consider robust to saturation and limiting 
signal-to-noise (red data points). Also shown on the figures 
are the derived global offset values (blue lines), and the sim¬ 
ple linear fits (green lines) to the medians (black squares 
with errorbars) for both the magnitude (left) and colour 
comparisons (right). The dotted lines indicate the quartiles 
of the data. 

We conclude that the absolute zero-point calibration is 
robust across the board to ±0.02 mag within the magnitude 
and colour ranges indicated (red data points). However, it 
is extremely important to recognise that the majority of our 
galaxies lie significantly outside the flux and colour ranges 
which we are examining here. As we shall discuss in Sec¬ 
tion this can cause significant and intractable issues. To 
quantify the potential for zero-point drift between the cal¬ 
ibration regime and operating regime we show in Table 
possible zero-point offsets one might derive at the typical 
flux and the typical colour of the GAMA sample using ei¬ 


ther (a) a simple offset (Fig.f^V 17 blue line), (b) a linear 


fit with magnitude (the linear fit shown as a green line in 
Fig. V 17 left panels) and, (c) a linear fit with colour 
(the linear fit shown as a green line on Fig. 16 V 17 right 
panels). Any one of these relations, or some combination of, 
could be valid and hence the range reflects the uncertainty in 
the absolute zero-point calculations for our filters. We elect 
not to correct our data using any of these zero-points but 
instead incorporate the possible systematic zero-point error 
(indicated in the final column) into our analysis. 


3.3 u — Kg photometry errors 


Gritical to any SED fitting algorithm will be the derivation 
of robust errors for each of our galaxies in each band. Here 
we derive the errors from consideration of: the zero-point 
error (crzp), the random sky error (crskyRan), the systematic 
sky error (crskySys), and the object shot noise (crshot)- The 
first of these is quoted in Table the other three can be 
given by: 


C^SkyRan = V N’pixCTSky 
^Sky 


C^SkySys = A^Pix 




C^Shot — 


Ap er 
-^Obj 

7 


(5) 

(6) 

(7) 


Where crsky is the sky noise give in Table (Gol.4), Vpix 
is the number of pixels in the object aperture (given by 
'^T^kron^imAGE Bimage in terms of Source Extractor out¬ 
put parameters), and AAper is the number of pixels used in 
the aperture in which the local background was measured 
(i.e., ttBkr^on (^iMAGE ± 32)(Bimage ± 32) — Vpix) and 7 is 
the gain. Of these only the gain is uncertain as during the 
stacking and renormalising of the data the gain is modified 
from its original value by varying amounts (see for example 
the distribution of multipliers in Eig.|^. However, as the vast 
majority of our galaxies are relatively low signal-to-noise de¬ 
tections the sky errors swamp the object shot noise errors 
and hence we elect to omit the object shot noise component 
in our final error analysis. 


3.4 Comparison to earlier GAMA photometry 


Einally we compare our revised SDSS±VIKING photome¬ 
try to our earlier SDSS±UKIDSS photometry in Eig. 

In this implementation of IOTA the only difference in the 
ugriz bands is the move from a global background estima¬ 
tion (fixed value across the background subtracted frame) 
to a local background estimation. The impact appears min¬ 
imal with zero-point offsets less than ±0.015. However in 
the Y JHKg bands we notice significant offsets between the 
UKIDSS and VIKING flux measurements. The reasons for 
this are subtle and while they have not been exhaustively 
pursued we believe are most likely due a hidden linearity is¬ 
sue in the UKIDSS pipeline. In Eig.[^we show our flux mea¬ 
surements from our UKIDSS data for 420k SDSS selected 
stars for which we have 2MASS photometry. The agreement 
is once again good, however in all cases there are significant 
gradients in the data and significantly stronger than those 
we saw in the VIKING data (r.f. Eig. 17). Extrapolating the 
linear fits to the flux and colour regions where the majority 
of our galaxies lie we infer the level of offsets seen in Eig.p^ 
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Figure 16. Comparison between SDSS PSF photometry versus SDSS (IOTA) for 420k stars. Highlighted in red are data deemed to lie 
in the flux and colour regions for which comparisons can be made. 


Table 5. Zero-Point uncertainties in each band. 


Band 

GAMA Median 

GAMA Median 

Potential zero-point (ZP) offsets 

ZP unc. 

Adopted 


flux limit (mag) 

(J — Ks) colour (mag) 

Absolute 

linear with mag 

linear with colour 

Mean ± Std. 

ZP error 

u 

21.48 

0.37 

+0.011 

+0.018 

-0.020 

+0.003 ± 0.020 

0.02 

g 

20.30 

0.40 

+0.017 

+0.035 

+0.048 

+0.033 ± 0.016 

0.05 

r 

19.35 

0.41 

+0.018 

+0.032 

+0.032 

+0.027 ± 0.008 

0.03 

i 

18.88 

0.41 

+0.016 

+0.033 

+0.033 

+0.027 + 0.008 

0.03 

z 

18.60 

0.41 

+0.019 

+0.036 

+0.035 

+0.030 ± 0.008 

0.03 

Z 

18.61 

0.41 

+0.006 

+0.039 

+0.076 

+0.040 ± 0.035 

0.08 

Y 

18.38 

0.42 

NA 

NA 

NA 

NA 

0.10 

J 

18.16 

0.43 

+0.004 

+0.028 

+0.037 

+0.023 + 0.017 

0.04 

H 

17.84 

0.44 

-0.015 

-0.050 

+0.010 

-0.018 + 0.030 

0.05 

K 

17.69 

0.44 

-0.010 

+0.026 

+0.018 

+0.011 + 0.019 

0.03 
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Figure 17. Comparison between 2MASS stellar photometry versus VIKING (IOTA) for 420k stars. Highlighted in red are data deemed 
to he in the flux and colour regions for which comparisons can be made. 


The implication is that there may be a linearity issue with 
the UKIDSS calibration. Note that as Hill et al. (2010) has 
shown our in-house UKIDSS photometry agrees extremely 
well with that provided from the UKIDSS archive. We do 
not explore this issue further but, as a number of earlier 
GAMA papers are based on UKIDSS photometry, we in¬ 
clude the UKIDSS SWarps in the public release, while cau¬ 
tioning against their use. 


3.5 Optical motivated far-IR photometry 

To derive our far-IR photometry for every GAMA target we 
implement an optically-motivated approach (also referred to 
as forced-photometry). This technique closely follows the ap¬ 
proach developed by Bourne et al. (2012) for the Herschel- 
ATLAS team and which has been used to obtain SPIRE 
photometry at the location of known optical sources. The 
method adopts as its starting point the r-band apertures 
determined from our optically motivated source finding de¬ 
scribed earlier and uses the following parameter set for the 
apertures: right ascension, declination, major axis, minor 
axis and position angle. For each far-IR band the aperture 
defined by these parameters is combined with the appropri¬ 
ate PSF for each of the five bands (supplied by the Herschel- 
ATLAS team). The resulting 2D distribution therefore con¬ 
sists of a flat pedestal (within the originally defined aperture 


region) with edges which decline as if from the peak of the 
normal PSF. This soft-edge aperture can be imagined as a 
2D mesa-like distribution function which can now be con¬ 
volved with the data at the appropriate astrometric location. 
In the event of two mesas overlapping the flux is shared ac¬ 
cording to the ratio of the respective mesa functions at that 
pixel location, i.e., the flux is distributed using PSF and 
aperture information only and ignoring the intensity of the 
central pixel. Enhancements of this methodology are under 
development (Wright et al. 2015) and will include consider¬ 
ation of the central peak intensity along with the inclusion 
of interlopers (i.e., high-z targets), and iterations. The re¬ 
covered fluxes at this moment contain flux from the object, 
plus from any low-level contaminating background objects. 
To assess the level of contamination we made measurements 
in ~ 30, 000 apertures of comparable sizes to our object dis¬ 
tribution which were allocated to regions where no known 
5cr Herschel-ATLAS detection exists nor any GAMA object. 
The mean background level in these regions was found to be 
zero in all SPIRE bands, as expected since the maps are 
made to have zero mean flux and residual large scale emis¬ 
sion has been removed via the nebuliser step. In the PAGS 
data small background values were found as shown in Ta¬ 
ble. To correct for this effect the final step is to subtract 
the background values for the PAGS data using the effeetive 
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Figure 18. Comparison between GAMA ApMatchedv03 and GAMA ApMatchedvOG (i.e., PDR) ugrizZYJHKg galaxy photometry. 
Filter transformations as indicated in the text. 


Table 6. Sky background levels per pixel as derived from blank 
apertures. 


Sky background (Jy) 


Band 

G09 

G12 

G15 

100p,m 

0.0002669 

0.0001491 

0.0002739 

160/rm 

0.0002044 

0.0002436 

0.0003784 


aperture pixel number and factoring in shared pixels where 
apertures overlap. 

Figure compares our aperture-matched photome¬ 
try against the Herschel-ATLAS 5cr catalogue produced by 
Smith et al. (2012). In general the data agree reasonably 
well with offsets at the levels of 0.05dex (12%). 

There have been changes in the PACS calibration and 
map-making algorithms between the generation of the values 
in the Smith et al. 2012, Rigby et al. 2011 catalogues and this 
work, and these changes have been substantial. Offsets at 
this level are consistent with these changes. For the SPIRE 
data, there have been no significant changes to the calibra- 
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Figure 20. A comparison between the preliminary Herschel- 
ATLAS 5(7 catalogue and the optically motivated catalogue de¬ 
rived here in the SPIRE 250p,m band. It should be noted that the 
zero-points calibration and entire reduction process have evolved 
between these catalogues. Note units are in Jansky as is standard 
in the far-IR. 



tion or map-making process, however there are a number of 
potential issues with both the catalogues being compared in 
in Fig |20[ Firstly, the H-ATLAS catalogue is with a prelimi- 
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Figure 19. Comparison between 2MASS stellar photometry versus UKIDSS (IOTA) for 420k stars. 


nary version of the H-ATLAS release catalogue (Valiante et 
ah in prep) which does not include aperture photometry for 
resolved sources, explaining some of the scatter at the bright 
end. Secondly, the largest optical sources, which correspond 
to the brightest H-ATLAS sources are often shredded by 
SExtractor which leads to inappropriately small apertures 
being used for the forced photometry, and this may lead 
to the offset between the catalogues at the bright end and 
contribute to the scatter. 

In Wright et al. (2015) we will compare our updated 
LAMBDAR photometry with the final released version of 
the H-ATLAS catalogue (Valiante et al. 2015) when all 
fluxes will be drawn from the same data pipelines and im¬ 
ages. 

3.6 Table matching the UV, optical/near-IR, 
mid-IR and far-IR catalogues 

At this stage we have a number of distinct catalogues. 

GalexMainv02: This contains measurements of the FUV 
and NUV fluxes which have been assembled through the 
use of r-band priors combined with curve-of-growth analysis 
and is described in Liske et al. (2015). We adopt the BEST 
photometry values. In brief the BEST photometry is that re¬ 
turned by the curve-of-growth method with automatic edge- 
detection when the NUV semi-major axis is greater than 20 
arcsec or when the GAMA object does not have an unam¬ 


biguous counterpart. In other cases the BEST photometry 
is that derived from the standard pipeline matched to the 
GAMA target catalogue. 

ApMatchedCatvOG: This contains the u to Kg band pho¬ 
tometry as described in detail in Section [3.1| 

WISEPhotometryv02: as described in Gluver et 
al. (2014) which outlines the detailed construction of the 
WISE photometry with two exceptions. Firstly, for GAMA 
galaxies not resolved by WISE, standard aperture photom¬ 
etry (as provided by the AllWISE Data Release) is used 
instead of the profile-fit photometry (wpro). This is due to 
the sensitivity of WISE when observing extended, but unre¬ 
solved sources, resulting in loss of flux in wpro values com¬ 
pared to standard aperture values (see Gluver et al. 2015 
for details). Secondly, the photometry has been updated to 
reflect the AllWISE catalogue values. Note that this version 
of the catalogue also includes the correction to the updated 
W4 filter described in Brown, Jarrett V Gluver (2014). 

HAtlasPhotomCatvOl: This contains the far-IR mea¬ 
surements as described in Section [375| based on optically mo¬ 
tivated aperture-matched measurements incorporating con¬ 
tamination corrections. 

We use TOPCAT to combine these catalogues by matching 
on GAMA GATAIDS (i.e., exact name matching), the FUV 
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to Ks-hand data are then corrected for Galactic extinction 
using the E{B — V) values provided by Schlegel, Finkbeiner 
& Davis (1998; GalacticExtinctionv02) and the coefficients 
listed in Liske et ah (2015). 

The combined catalogue is then converted from a mix¬ 
ture of AB mags and Janskys to Janskys throughout, with 
dummy values included when the object has not been sur¬ 
veyed in that particular band. Goverage maps may be re¬ 
covered from the catalogue using the dummy values alone. 


4 ROBUSTNESS CHECKS OF THE PDR 


As the GAMA PDR is constructed from a variety of distinct 
catalogues and pathways it is important to assess its robust¬ 
ness, accuracy, and outlier rate. In earlier figures (Figsto 
19) we showed direct comparisons of the magnitude differ¬ 
ence between two datasets. These are good for identifying 
zero-point (i.e., systematic) offsets, but not particularly use¬ 
ful in establishing which of the two datasets is the more ro¬ 
bust. Here we examine the more informative “colour”-plots. 
The implicit assumption is that a colour distribution arises 
from a combination of the intrinsic colour spread of a galaxy 
population, convolved with the measurement error in the 
contributing filters. A comparison of colour-plots between 
two surveys, for the same sample, can provide two impor¬ 
tant statistics: the width of the distribution, and an outlier 
rate. The “better” quality data is the dataset with the nar¬ 
rowest colour range and the lowest outlier rate (assuming 
the zero-points are consistent). The colour-plot test is opti¬ 
mal when the intrinsic colour spread is sub-dominant, hence 
should be made using adjacent filters. In some bands, e.g., 
(NUV — u) the intrinsic colour range is known to be broad 
(e.g., Robotham & Driver 2011), and hence the test less 
conclusive. Fig. shows the full set of colour distributions 
for the GAMA PDR (black histograms), the data have been 
Galactic extinction corrected but not /c-corrected, and this 
is chosen to minimise the modeling dependence, particularly 
given the wavelength range sampled and uncertainty in k- 
correcting certain regimes (e.g., mid-IR). 

Also shown on Fig. is the breadth of the colour dis¬ 
tribution derived from the 80%-ile range (horizontal red line 
and red text) and the median colour value (vertical red line). 
We derive an outlier rate (indicated by “Out” as a per¬ 
centage on the figure), this reports the percentage of galax¬ 
ies which lie more then 0.5mag outside the 80 percentile 
range. The rationale is that the 80 percentile distribution 
will generally capture the intrinsic+Zc-correction spread, and 
a catastrophic magnitude measurement would then be one 
which lies more than 0.5mag outside of this range. One can 
see that the colour distributions are particularly broad in the 
UV bands (as one expects given the range from star-forming 
systems to inert systems with varying dust attenuation), and 
in the far-IR bands (as one expects given the range of dust 
masses and dust temperatures). The red optical and near-IR 
bands are the narrowest (as expected given the flatness of 
SEDs at these wavelengths). The outlier rates are generally 
highest for the poorest resolution and lowest signal-to-noise 
bands (i.e., NUV, FUV, u and W2 onwards) with outlier 
rates varying from 12.4 percent to 0.5 percent. Our ultimate 
objective within the GAMA survey is to achieve outlier rates 
below 2 percent in all bands. With 10 colour distributions at 


or below this level this implies we have reached this criterion 
for 11 bands {g to Wl). 

The facility cross-over colours, {NUV—u), (z—Z), {K— 
VFl), (1F4—lOO/xm), and (IGO/xm- 250/im), are of particular 
interest as this is where mis-matches between objects might 
lead to broader distributions and higher outlier rates, and 
there is some indication that outlier rates do rise at these 
boundary points, e.g., the {NUV — u) and {z — Z) bands 
and (K — Wl) bands. In the case of the former this may 
simply reflect intrinsic+/c-correction spread. 

No dataset sampling a comparable wavelength range 
currently exists. However, we can compare in the optical 
and near-IR to the SDSS archive, our previous catalogue 
(based on SDSS and UKIDSS) and into the UV and mid- 
IR with the low-z templates given in Brown et al. (2014). 
These are shown where data exists as blue (GAMA Ap- 
MatchedGatvOS; SDSS+UKIDSS), green (SDSS DR7 Mod- 
elMags), and purple histograms (Brown et ah). Note that 
as the SDSS and GAMA data are essentially derived from 
the same base optical data and it is the photometric mea¬ 
surement method which is being tested here. In future we 
will be able to compare to KiDS and Subaru HSG datasets. 
In comparison to our previous GAMA catalogue we can see 
that PDR represents an improvement (lower breadths and 
lower outlier rates) in all bands. In particular the near-IR 
bands are significantly improved with the colour spread now 
at least two times narrower. This reflects the greater depth 
of the VIKING data over the UKIDSS LAS data (see Ta¬ 
ble]^. In comparison to SDSS DR7 ModelMags we can see 
that GAMA PDR appears to do marginally better in the 
(u — g) and {g — r) colours, but marginally poorer in the 
(r — i) and {i — z) bands, in all cases by modest amounts. 

Einally, in comparison to the colour distributions de¬ 
rived from the Brown et al. (2014) templates there are two 
important caveats. Eirstly the Brown data makes no attempt 
to be statistically representative but rather provides an in¬ 
dication of the range of SEDs seen in the nearby population 
for a relatively ad hoc sample. Secondly the GAMA PDR 
has a median redshift of z = 0.24 and in some bands the k- 
correction will dominate over the intrinsic distribution. This 
is apparent in particular in the (NUV—u), (u—g) and (g—r) 
bands where the 4000A break is redshifted through. This re¬ 
sults in significantly broader colours in the observed GAMA 
PDR colour distributions not seen in the rest-frame tem¬ 
plates. Again this is understandable. More puzzling is the 
converse where the (W2 — Wd) and (IU3 — WA) distribu¬ 
tions which are clearly broader in the Brown et al data. This 
may reflect the incompleteness within the GAMA PDR in 
these bands, with the bluest objects perhaps being detected 
in W2 but not in W3, and hence not represented on these 
plots (see Gluver et al. 2014 for full discussion on the WISE 
completeness). Similarly for {W3 — W4). Eor example the 
Brown et ah, sample includes both Elliptical systems and 
very low luminosity blue dwarf systems (e.g., Mrk 331, H 
Zw 96, Mrk 1490, and UM 461) neither of which would be 
likely detected by WISE at z >> O.OI. The obvious solution 
is to derive “forced-photometry” for the full GAMA input 
catalogue across all bands. 

At this point we believe we have established that 
GAMA PDR is matching SDSS DR7 ModelMags, a sig¬ 
nificant improvement over previous GAMA work based on 
SDSS+UKIDSS LAS, but there remains some concerns re- 
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garding higher than desired outlier rates in the lower signal- 
to-noise and poorer resolution bands, and the need for a 
measurement at the location of every GAMA galaxy regard¬ 
less of whether there is obvious flux or not (i.e., forced pho¬ 
tometry). Fixing these problems is non-trivial and requires 
dedicated panchromatic software, which is currently nearing 
completion and will be presented in Wright et ah (2015). 


4.1 Composite SEDs 

Using a 35,712 core machine available at the Pawsey Su¬ 
percomputing Centre Facility (MAGNUS) we have now run 
the MAGPHYS spectral energy distribution fitting code 
(de Cunha, Chariot & Elbaz 2008), over the full equatorial 
GAMA sample with redshifts, i.e., 197k galaxies (using the 
Bruzual & Chariot 2003 spectral synthesis model). MAG¬ 
PHYS takes as its input, flux measurements in each band, 
associated errors, and the hlter bandpasses, and returns the 
attenuated and unattenuated SED models from EUV to far- 
IR, along with a number of physical measurements, e.g., 
stellar mass, star-format ion rate, dust mass, dust opacity, 
dust temperatures (birth-cloud and ISM) etc (for more de¬ 
tails please see da Cunha et al. 2008). Here we look to use 
MAGPHYS to provide a simple spectral energy represen¬ 
tation for each of our galaxies which also has the effect of 
hlling in the gaps where coverage in a particular band does 
not exist or no detection is measured. On a single processor 
MAGPHYS will typically take lOmins to run for a single 
galaxy (i.e., 4yrs for our full sample), but using MAGNUS 
the entire sample can be processed in less than 24hrs. 

Eig.[^ shows each of the 21 x 197,491 datapoints plot¬ 
ted at the rest wavelength versus the number of a devia¬ 
tions the GAMA PDR magnitude is away from the derived 
MAGPHYS magnitude. The figure highlights that generally 
MAGPHYS appears to be hnding consistent fits across all 
bands with only the W4 showing some indication of a fun¬ 
damental inconsistency between the data and the models. 
This has now been tracked down to the change in the W4 
hlter transmission curve with our MAGPHYS run still using 
the old throughput curve while the WISE PDR data uses 
the revised W4 transmission curve (see Brown, Jarrett & 
Cluver 2014 for full details). Euture runs of MAGPHYS will 
use the updated curve. The distribution of the data points in 
a deviations (abscissa) suggest that the measurement vari¬ 
ations are consistent with the errors quoted. The bleed of 
the far-IR data to the lower part of the hgure, is most likely 
due to contamination by high-z systems (as expected). In 
the wavelength range, within each hlter the near-IR data 
shows the most huctuations. These are likely to rehect re¬ 
current features in the MAGPHYS models shifting through 
the various bands and suggests some uncertainty in the pre¬ 
cise modelling of the TP-AGB region as noted by numerous 
groups, e.g., Maraston et al. (2006). 

Eig. shows the MAGPHYS SED model hts sampled 
by our z < 0.06 morphologically classihed sample, sepa¬ 
rated into E/SOs (red), Sabcs (green) or Sd/Irrs (blue). 
See Moffett et al. (2015) for details on the sample and 
morphological classihcation process. In order to construct 
these plots the individual SEDs derived by MAGPHYS have 
been normalised to the same stellar mass, i.e., their SEDs 
have been scaled by their fractional stellar mass offset from 
10^°M* (using the stellar masses derived by MAGPHYS). 


The curves shown are the quantile distributions for 10, 25, 
45, 50, 55, 75 and 90 percent as a function of wavelength. 
Note that these do not represent individual MAGPHYS SED 
models, but are quantile ranges in narrow wavelength inter¬ 
vals which are then linked to create the SED quantiles - 
hence the SEDs show more variations than the models used 
in MAGPHYS if examined in detail. The reason for this 
representation is to avoid a specihc calibration wavelength 
and SEDs that, when calibrated into quantiles at one wave¬ 
length, cross at others. 

Eig. illustrates not only the wealth of data provided 
by the GAMA PDR but a number of physical phenomena. 
Eirstly the spread at any wavelength point represents the 
mass-to-light ratio at that wavelength. This can be seen to 
be narrowest in the 2 - 5/im range (as expected), indicating 
that this region is optimal for single band stellar-mass esti¬ 
mates. However, the constant gradient in this region implies 
that near-IR colours provide little further leverage to im¬ 
prove the stellar-mass estimation beyond single band mea¬ 
surements. Gonversely the smooth variation of SED gradi¬ 
ents in the optical from low to high stellar-mass ratios, im¬ 
ply that optimal stellar-mass estimation may arise from the 
combination of a single band near-IR measurement com¬ 
bined with an optical colour (see also discussion in Taylor 
et al. 2011). Eig. 23 highlights the known strong correla¬ 
tions between UV flux, far-IR emission and stellar mass-to- 
light ratio with all being amplified or suppressed in Sd/Irrs 
or E/SOs respectively. However, that all galaxies seem to 
contain some far-IR emission may be a manifestation of 
the MAGPHYS code tending to maximise dust content 
within the bounds as allowed by the far-IR errors. Guri- 
ously the Sabc (green) provide a very narrow range of pa¬ 
rameters, perhaps indicating a close coupling between star- 
formation, dust production, and the mass-to-light ratio — 
arguably indicative of well-balanced self-regulated disk for¬ 
mation/evolution. The greater spreads in the early (red) and 
later (blue) types are perhaps indicative of the progression 
through various stages of quenching (ramping down) and 
unstable disc formation (overshoot) respectively. In particu¬ 
lar Agius et al. (2015) found an unexpectedly high levels of 
dust in a significant (29%) population of the GAMA-E/SO 
galaxies consistent with a range of E-So SEDS. Note also 
the results presented on observed correlations between the 
star-format ion rate, specihc star-format ion rate in da Gunha 
et al. (2010) and Smith et al. (2012); see also interpretation 
in Hjorth et al. (2014). A detailed exploration of these phe¬ 
nomena are beyond the scope of this paper but the potential 
is clear particularly in conjunction with the existing group 
(Robotham et al. 2011) and large scale structure (Alpaslan 
et al. 2014) catalogues. 


4 . 1.1 Inspection of individual objects 

We explore individual SEDs for one hundred systems ran¬ 
domly selected (IDs 47500-47609). Approximately 5-10 per¬ 
cent are found to have one or more signihcant outlier(s) 
in the photometry but otherwise good MAGPHYS hts are 
found for all 100 systems. In approximately 50 percent of 
cases the far-IR photometry is essentially missing (due to 
the shallowness of the far-IR data), hence hux and/or red- 
shift cuts are advisable depending on the science investiga¬ 
tion to be conducted. Eig. shows four example galaxies 
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Figure 21. Each panel shows a histogram of a GAMA PDR colour histogram (black distribution) from two adjacent filters (as indicated, 
y-axis). The horizontal red bar indicates the 80-percentile range of the data with the value for GAMA PDR indicated in black in the top 
left corner. The outlier rate (indicated by Out) reports the percentage of the distribution which lies >0.5mag outside the 80-percentile 
range (indicated by the vertical dotted lines). Where possible we overlay data from the SDSS DR7 (ModelMag colours, green histograms) 
and ApMatchedGatvOS (which uses SDSS and UKIDSS data, blue histograms). All distributions are extinction corrected but otherwise 
as observed. Finally the shaded histogram (purple) shows data from the z = 0 sample of Brown et al. (2014). However, note that the 
Brown sample has an ad hoc selection and essentially represent k-band corrected data - hence explaining the apparent discrepancy in 
{NUV — n), {u — g) and {g — r), and {H — K). 
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Figure 22. Sigma deviations of the model photometry compared to the input photometry versus rest wavelength for all data points for 
all galaxies. For each colour the median and Icr dispersions is measured and shown (grey solid and dotted lines respectively). The plot 
shows that the majority of our data lies within the quoted error of the MAGPHYS fits with some outliers (as also noted in Fig. |21| ). The 
offset of W4 is of some concern as is the bleeding of the far-IR data towards negative deviation values. Both of these effects are believed 
to be understood but requires the LAMBDAR software (in development). If these are flat and offset it implies a systematic offset, if 
they show trends with wavelength it implies a progression of a feature with redshift. Trends are most apparent in the near-IR where the 
modelling of the TP-AGB population is still uncertain. 
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Figure 23. The panel shows the 10, 25 ,45 ,50 ,55 ,75 and 90 percentiles for the full GAMA PDR SEDs for the 2 ; < 0.06 morphologically 
classified sample presented in Moffett et al. (2015). The data are initially scaled to the same mass and then at each wavelength point 
the quantiles derived which then collectively trace out the quantiles over wavelength. The spread at each wavelength therefore directly 
reflects the spread in the mass-to-light ratio at that wavelength. 
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which include a nearby bright system (G47152, z = 0.082), 
a nearby faint system (G47157, z = 0.074), a higher red- 
shift crowded system (G47609, 2 : = 0.282), and a known far- 
IR lens system (G622892, 2 : = 0.300; Negrello et al. 2010; 
recently shown to exhibit a spectacular Einstein ring, the 
very high far-IR flux is evident). The panels on the left 
show the combined giH colour image from a combination 
of VIKING and SDSS data. Overlaid (green dotted lines) 
are apertures for the main object and nearby systems in 
our bright catalogue. The right panels show the 21-band 
measured photometry (green data and errorbars) in units 
of total energy output (XLx in units of hyo W) at the filter 
pivot-wavelength divided by (1 + z) (i.e., rest wavelength). 
The red and blue lines show the attenuated and unatten¬ 
uated SEDs from the preliminary MAGPHYS fits. Purple 
circles show the flux from the attenuated SED curve inte¬ 
grated within the filter bandpasses given in Eig.[^ The lower 
portion of the panel shows the residuals expressed as the ra¬ 
tio of the observed flux to the measured flux. Included in 
the error budget is a 10 percent flux component added in 
quadrature to mitigate small systematic zero-point offsets 
at facility boundaries. Gomparable plots for all 221k sys¬ 
tems with redshifts are provided via the GAMA T online 
cutout tool (http://gama-psi.icrar.org/) 

The four galaxies are each well sampled by the GAMA 
PDR data which collectively map out the two key peaks in 
the energy output due to starlight and dust reprocessing of 
starlight. The four systems also show varying degrees of dust 
attenuation with the observed fluxes requiring minimal, sig¬ 
nificant and extreme corrections to recover the unattenuated 
fluxes. In all cases the residuals are well behaved, given the 
errors. Exploring the GAMA PDR more generally, the SED 
data appear robust with catastrophic failures in one band 
typically at the levels indicated in Eig.[^ i.e., 10 percent in 
poorly-resolved bands to 1 percent in well-resolved bands. 
Obvious issues which arise following inspection of several 
hundred SEDs are: incorrect apertures, data artifacts (i.e., 
poor quality regions), nearby bright stars (diffraction spikes 
and blocking), crowding, and confusion. 


5 THE ENERGY OUTPUT OF THE 

UNIVERSE AT Z < 0.2 FROM FUV TO 
FAR-IR 

We conclude this paper with a brief look at the integrated 
energy output of the low redshift galaxy population, i.e., 
the Gosmic Spectral Energy Distribution (GSED), its re¬ 
cent evolution, and the implied integrated photon escape 
fraction. The GSED represents the energy output of a cos- 
mologically representative volume, in essence an inventory of 
the photons recently generated, as opposed to those passing 
through but formed earlier. It can be reported both pre- 
and post- attenuation by the dust content of the galaxy 
population, both are interesting. The pre-attenuated GSED 
informs us of the photons being created from (primarily) 
nucleosynthesis processes (in the current epoch), while the 
post-attenuated GSED informs us of the photons entering 
into the IGM. The sum of the two must equal (energy con¬ 
servation) , but the wavelength distribution will differ as dust 
re-processes the emergent photons from short (UV and op¬ 
tical) to long (mainly far-IR) wavelengths. The combina¬ 


tion of the two can be used to determine the integrated 
photon escape fraction. By integrated we imply over a rep¬ 
resentative galaxy population, and averaged over represen¬ 
tative viewing angles. Both of these factors are important 
and make the integrated photon escape fraction useful for 
converting observed EUV and NUV fluxes to robust star- 
formation rates. The work follows earlier measurements of 
the GSED reported in Driver et al. (2008) and Driver et 
al. (2012). However, the methodology here is very different 
and for the first time includes mid- and far-IR data in a fully 
consistent analysis. In our earlier studies we determined lu¬ 
minosity distributions in each band independently and then 
fitted across these values to determine the GSED. Here we 
stack the individual MAGPHYS SEDs fits derived earlier, 
as representative fitting functions (see Eig. 24). 

Potentially, as we have a MAGPHYS fit for every galaxy 
we could use them to derive fluxes in data gaps and use the 
full sample. However, given the critical importance of the 
far-IR dust constraint we elect to use just the common re¬ 
gion with full 21-band coverage (see Eigs. ^^to 2 ^. This 
combined region constitutes an area of 63 percent of the full 
area or 113 deg^ and contains I38k objects with secure red- 
shifts in the range 0.02 < 2 : < 0.5 (trimmed to exclude 
stars and high-z AGN). To explore any evolution of the 
GSED, we divide our sample into three redshift intervals: 
0.02 < z < 0.08, 0.08 < z < 0.14 and 0.14 < z < 0.2 
which correspond roughly to lookback times of 0.8 Gyr, 
1.5 Gyr and 2.25 Gyr respectively. The volumes sampled are: 
4.9x 10^, 2.1 X 10® & 4.6x 10® h^Q® Mpc® respectively (fac¬ 
toring in our reduced coverage). Within each redshift range 
we use the Zmax values reported in Taylor et al. (2011) to 
derive a weight as not all galaxies would be visible across the 
selected redshift range. Galaxies with Zmax values above the 
redshift range have weights set to unity, and values with Zmax 
below this range have weights set to zero. Otherwise weights 
are set to the inverse of the fraction of the volume sampled. 
A cap is placed ensuring no weight exceeds a value of 10, 
this ensures a single lone system just fortuitously within 
the redshift range cannot dominate the final outcome by 
being massively amplified. Within each redshift range we 
now simply sum the energy x weight (i.e., J^H5AL\) for 
the galaxies within our selection to arrive out the GSED for 
that volume. 

These raw derived GSEDs require one final correction 
to accommodate for the loss of lower luminosity systems in 
the higher redshift bins. To determine the correction fac¬ 
tor we repeat the summation but with a stellar mass cut 
imposed on all three samples (I0^®Mq), this is sufficiently 
low to be sampling the dominant contribution to the GSED 
but not so low as to suffer total incompleteness (i.e., that 
which is not corrected for by our weights). Eor each volume 
interval we obtain unrestricted to mass-restricted GSED ra¬ 
tios of: 1.68, 1.48 and 1.29 for the low, mid, and high red¬ 
shift samples respectively. If all three samples were complete 
this ratio would be constant, hence this changing ratio en¬ 
codes the loss of the lower luminosity systems in the higher 
redshift bins, and can therefore be used to provide an ap¬ 
propriate correction. This is achieved by scaling the final 
GSED curves by factors of 1.00, 1.14 (i.e., 1.68/1.48), and 
1.30 (i.e., 1.68/1.29) respectively. In effect we are using the 
GSED shape from the unrestricted samples but normalis¬ 
ing using the restricted samples and this is analogous to the 
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Figure 24. Four example galaxies from the GAMA PDR showing, in each case, the Hig image with the apertures overlaid, and the 
21-band photometry with the attenuated (red) and unattenuated (blue) MAGPHYS fits to the data. Data shown in green represent 
the measured photometry plotted at rest wavelength while the purple circles show the integrated fiux measured from the attenuated 
MAGPHYS fit through the appropriate filter. The lower portion of the main panels shows the residuals, i.e., the ratio of the observed 
fiux to the model fiux. Errors include an arbitrary 10 percent error added in quadrature intended to incorporate some allowance for 
zero-point offsets between facilities. Note that the lowest image is a well known lens system reported by Negrello et al. (2010). All images 
are derived from the GAMA T online cutout tool http://gama-psi.icrar.org/ 
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normalisations typically used in estimating luminosity and 
mass functions. This implicitly assumes the following: that 
the low redshift sample is itself complete (hence requiring 
no scaling), and that the ratio of energy emerging from sys¬ 
tems above and below is approximately constant. 

The first of these is relatively secure: GAMA is a deep sur¬ 
vey and at 2 : ~ 0.08 is mass complete to lO^M© (see Lange et 
al. 2015), below which there is very little contribution to the 
luminosity density (see Driver 1999), or stellar mass density 
(see Moffett et al. 2015). In the second case we understand 
low-mass systems are preferentially star-forming and may 
have evolved more rapidly over recent times compared to the 
more massive systems. However, as the correction factors are 
relatively modest (14 percent and 30 percent) the shape and 
renormalisation is unlikely to be dramatically changed, but 
we acknowledge may be biased low. This can only be quan- 
tihed through deeper studies (see for exampled the planned 
WAVES survey; Driver et al. 2015). 

Fig.|25|shows the resulting energy outputs for the renor¬ 
malised unattenuated CSED (upper) and the renormalised 
attenuated CSED (lower) with the redshift ranges repre¬ 
sented by colour as indicated. Also shown is our earlier 
estimate derived from GAMA via luminosity function fit¬ 
ting (orange line). It is worth re-iterating that the Driver 
et al. (2012) CSED measurement is based on luminosity 
density measurements from FUV to Ks combined with 
an adopted integrated photon escape fraction (Driver et 
al. 2008) to infer the mid and far-IR portion, i.e., the mid 
and far-IR from Driver et al. (2012) is a prediction, and 
hence shown as a dotted line beyond 2.1/im. Also shown is 
the prediction from semi-analytic modelling by Somerville 
et al. (2012). Both curves, corrected to Ho = 70 km/s/Mpc, 
follow the new low-z CSED very well in the optical and 
start to diverge in the mid and far-IR bands where previous 
empirical data have been lacking. The Somerville curve, in 
particular traces the low-redshift bin extremely well with the 
previous Driver data significantly under-predicting the far- 
IR emission. The CSEDs presented here represent a major 
advance constituting the first consistent measurement of the 
post- attenuated CSED from a single sample spanning from 
the FUV to far-IR. Hence, while cosmic (sample) variance 
(CV) may scale the respective CSEDs in overall density, it 
will not modify the shape of the distribution (unless there 
are extreme hidden clustering factors). Using the formula 
in Driver & Robotham (2010, see also the online calculator 
http://cosmocalc.icrar.org/) we derive the CV, for our three 
redshift ranges to be: 18 percent, 12 percent and 10 percent 
respectively, with increasing redshift and note that despite 
the extent of the GAMA PDR the dominant error remains 
the CV. 

To first order the CSEDs are therefore all consistent 
with each other and the previous values. The CSEDs fol¬ 
low the expected progression towards higher energy output 
towards higher redshifts, despite the potential uncertainty 
in overall normalisation from CV. In particular the far-IR 
increases noticeably faster than the optical. This result is 
independent of cosmic variance and noted previously within 
the far-IR community, see for example Dunne et al. (2011) 
who infer a significant increase in dust mass (towards high 
redshift) over the past 5 Gyr. In detail we can report that 
the Universe is in energy decline, having dropped from a to¬ 
tal energy production of (2.5 ± 0.2) x 10^^ hro W Mpc“^ at 


2.25 Gyr ago, to (2.2 ± 0.2) x 10^^ hro W Mpc”^ 1.5 Gyr 
ago and (1.5 ±0.3) x 10^^ hro W Mpc“^ 0.75 Gyr ago. This 
decline is significant despite the CV uncertainty, and in line 
with our understanding of the evolution in the cosmic star- 
formation history (e.g., Hopkins V Beacom 2006), which 
shows a decline of a factor of approximately xl.5 over this 
time-frame. The lowest redshift bin value is also consistent 
with the (1.26 ±0.09) x 10^^ hro W Mpc“^ for z < 0.1 as re¬ 
ported in Driver et al. (2012) and the orange and blue curves 
on Fig.j^do show consistency over the optical regime. The 
majority of the energy difference is derived from the far- 
IR where previous data was lacking and the CSED in this 
region estimated. As an aside it is also worth noting that 
the near-IR data show some slight disagreement, and this 
is consistent with our finding that VIKING fluxes are typ¬ 
ically brighter than UKIDSS LAS fluxes (see Section 3.4). 
Finally, shown on Fig.j^are the SPIRE data points derived 
by the Herschel-ATLAS team from the initial data release 
which sampled all of GOO, G15 and half of Cl2 (Bourne 
et al. 2012). These data are low compared to our new es¬ 
timate, however the SPIRE calibration and data reduction 
have evolved quite substantially since these data-points were 
derived. 

Integrating these distributions reveals some interesting 
numbers. From the unattenuated CSED we hnd that 50 per¬ 
cent of the energy production in the nearby Universe is at 
wavelengths in the range 0.01—0.64/im (i.e., UV/optical). In 
the post-attenuated “observed” CSED, 50 percent of the en¬ 
ergy emerges at wavelengths shorter than 1.7/im. Splitting 
at lO/xm we hnd that 65 percent of the energy produced 
(via stellar emission) is released into the IGM at shorter 
wavelengths and 35 percent at longer wavelengths. Despite 
dust contributing a very small proportion of a galaxy’s mass 
(typically < 1 percent; Driver et al. 2008), its impact on the 
energy output is dramatic with signihcant potential conse¬ 
quences for optically based hux and size measurements (see 
for example Pastrav et al. 2013). 

Dividing the pre- and post-attenuated CSEDs yields the 
integrated photon escape fraction (IPEF). This is a par¬ 
ticularly interesting distribution as it encapsulates the im¬ 
pact of dust in a simple and general way. Fig. |26| shows the 
IPEF in the three redshift bins. Also shown (purple dia¬ 
monds) is the IPEF derived in Driver et al. (2008) for the 
Millennium Galaxy Catalogue (Liske et al. 2003). Clearly 
apparent is a trend towards lower photon escape fractions 
towards modestly higher redshift, demonstrating the impact 
of dust evolution in our perception of even low redshift sys¬ 
tems {z < 0.2). It is important to note that this result is 
resilient to CV as we are comparing the ratio of the pre- and 
post- attenuated CSEDs and hence the normalisation can¬ 
cels out. The variation appears smooth with redshift and 
significant, with the escape fraction in the FUV changing 
from 18 percent in the higher redshift bin to 23 percent in 
the lower z bin. The corresponding change in the NUV is 
27 percent to 34 percent respectively. If this trend continues 
the implication is that at even intermediate redshifts the UV 
photon escape fraction may be significantly lower than the 
values typically adopted when deriving star-format ion rates 
from UV fluxes. 

The above work should be considered as preliminary, 
but indicative of the potential of the GAMA PDR to ex¬ 
plore the energy outputs of galaxies and galaxy population 
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Figure 26. The photon escape fraction integrated over all view¬ 
ing angles and derived from Fig. |25| Also shown is the escape 
fraction reported by Driver et al. (2008) from the Millennium 
Galaxy Catalogue within the redshift range 0.0 < 2 ; < 0.18. The 
data show a clear progression towards lower escape fractions as 
redshift increases. 


at low redshift. The analysis at present also includes a num¬ 
ber of important caveats which we are looking to address 
in the near future. The first is out-standing issues related 
to the GAMA photometry: aperture robustness, the need 
for a uniform forced photometry across all bands, and im¬ 
proved management of the variable spatial resolutions and 
signal-to-noise limits. The second is whether MAGPHYS 
provides truly unbiased fits, particularly in the far-IR where 
our data quality is lowest and where, as the errors grow, 
over-fitting of the far-IR fluxes cannot be ruled out. Veri¬ 
fying that MAGPHYS is unbiased will be important as our 
data quality undoubtedly declines with redshift. One should 
also bear in mind that the interpretations above are very 
much based on the assumptions embedded in the MAG¬ 
PHYS code, in due course it will be important to explore 
a range of assumptions and to undertake critical compar¬ 
isons against fully radiative-transfer codes for well-resolved 
systems. The third and arguably most fundamental issue 
relates back to the integrity of the GAMA input catalogue 
and in particular its reliance on the fairly shallow SDSS 
imaging and the potential for missing extended low surface 
brightness systems. The GAMA regions are currently being 
surveyed by the VST KiDS team and will also be surveyed 
as part of the HSG Wide survey, both datasets can be used 
to improve our input catalogue particularly for extended low 
surface brightness systems. 


6 SUMMARY 

We have brought together a number of diverse datasets 
from three space missions (GALEX, WISE and Herschel), 


and two ground-based facilities (SDSS and VISTA) to pro¬ 
duce the GAMA Panchromatic Data Release. The individual 
dataframes have been astrometrically and photometrically 
matched, and then SWarped into single images for each of 
the three equatorial GAMA regions covering 80 deg^ each 
(i.e., slightly larger than the nominal 60 deg^ region covered 
by the spectroscopy). Weight maps are provided indicating 
the number of frames contributing to each region and hies 
containing the names of the individual frames used. In the 
u — Ks bands we provide both native resolution data and 
data degraded to a common 2" full-width half maximum. 
The SWarped images along with a cutout tool (GAMA ^) 
for extracting subregions are available at: 

http://gama-psi.icrar.org/ 

Note that the Herschel data is currently proprietary but 
will be made available following the Herschel-ATLAS hnal 
data release. 

GAMA ^ also provides additional functionality to cre¬ 
ate colour images on the hy from the combination of any 
of the 21-bands, individual hts downloads, object and aper¬ 
ture overlays and is maintained by the International Gentre 
for Radio Astronomy Research (IGRAR) Data Intensive As¬ 
tronomy (DIA) unit. Queries and comments on GAMA ^ 
should be sent to simon.driver@uwa.edu.au. 

We describe the construction of the u — Ks aperture- 
matched photometry following the method described in Hill 
et al. (2011) and compare to our previous measurements. In 
addition we compare our near-IR photometry for 400k stars 
extracted from 2MASS and conhrm that our zero-points are 
robust to within a tenth of a magnitude. 

The VIKING-2MASS and UKIDSS-2MASS zero-points 
are shown to be robust however VIKING-UKIDSS photome¬ 
try do not agree. This apparent tautological inconsistency is 
most likely explained by a linearity issue as the VIKING and 
UKIDSS data are compared at a signihcantly fainter mag¬ 
nitude range. Some evidence is seen for a linearity issue in 
the UKIDSS-2MASS comparison, which when extrapolated 
to fainter magnitudes does appear to explain the offsets seen 
between VIKING and UKIDSS. 

We also describe our method for deriving far-IR mea¬ 
surements using optical priors where we use our r-band aper¬ 
tures convolved with the appropriate instrument/hlter PSE 
to measure the far-IR hux at the location of every object. 

We then combine our u — Ks and 100 — 500/im huxes 
with GALEX and WISE data derived by the MPIK and 
UGT/UWG groups led by RJT and MG using exact name 
ID matching. This results in a hnal catalogue of over 221k 
objects from the GOO, G12, G15 and G23 regions, 63 percent 
of which have complete coverage in all 21 bands. 

Einally as a demonstration of this dataset we derive 
the total hux originating and emanating from various vol¬ 
umes from EUV to far-IR at redshifts indicative of 0.8, 1.5 
and 2.25 Gyr lookback time. We see evidence for evolution 
over this period consistent with the decline of star-formation 
traced by the cosmic star-formation history and consistent 
with the reported evolution of the far-IR luminosity func¬ 
tions. 

In future papers we will further improve our hux extrac¬ 
tion method using optically motivated priors in ah bands, 
leading to more consistent errors, and explore the physical 
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Figure 25. The energy originating (i.e., unattenuated, top), and emanating (i.e., attenuated following dust reprocessing, lower) at 
intervals equivalent to 0.75, 1.5 and 2.25 Gyr lookback time. The data are normalised to the energy output per Mpc^ for Ho = 
70km/s/Mpc. The data show clear trends in the evolution of the total energy output over this timeline. 


properties derived from SED fitting codes for various popu¬ 
lations and sub-populations. 
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APPENDIX A: GAMA PDR COVERAGE 

Eigs.[^to[^shows the GAMA PDR coverage in each of the 21 
bands, all data is available for download via: http://gama- 
psi.icrar.org/ 
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Figure 1. Background uniformity and coverage in the G09 region for GALEX, SDSS, VIKING, WISE and Herschel data. The black box 
denotes the GAMA spectroscopic survey region and the blue box the region from which the background statistics were derived. We use 
the MOGRIFY package to display the low-res data frames very close to the sky level. Astronomical objects will not be visible, however 
the frames highlight the coverage, missing regions, and the integrity of the large scale sky structure. 
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Figure 2. Background uniformity and coverage in the G12 region for GALEX, SDSS, VIKING, WISE and Herschel data. 


© 2010 RAS, MNRAS OOO.PHOOI 









































































GAMA Panchromatic Data Release 33 


z s-rfijaiTV 


tii.vn-,az 





:i'i 'j ir.:! h 


IK 





Figure 3. Background uniformity and coverage in the G15 region for GALEX, SDSS, VIKING, WISE and Herschel data. 
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Figure 4. Background uniformity and coverage in the G23 region for GALEX, VIKING, WISE and Herschel data. 
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